TOOL 

- ■  -  -  :■::   --:  ■>- .  .- :  -A  b J94S--600S 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey ,  California 


THESIS 


AN  ADVANCED  STUDY  OF  NATURAL 
ir'J-.ERSION  COOLING  OF  A  3  x 
SIMULATED  COMPONENTS  IN  AN 
FILLED  WITH  DIELECTRIC 

CONVECTION 
3  ARRAY  OF 
ENCLOSURE 
LIQUID 

by 

Terry  J.  Benedict 

June  1988 

Co- 

-advisors : 

Matthew  D.  Kelleher 
Yogendra  Joshi 

Approved  for  public  release;  distribution  is  unlimited 


12^8696 


UlNUijrtODXr   ±L.U 


SECURITY  CLASS. -CA'  ON 


n; 


■^■S    =AGc 


REPORT  DOCUMENTATION  PAGE 


la    REPORT  SECURITY  CuASSiFiCATiON 
UNCLASSIFIED 


lb    RESTRICTIVE   MARKINGS 


2a    SECURITY  CLASSIFICATION  AUTHORITY 


2b    DECLASSIFICATION 'DOWNGRADING   SCHEDolE 


3     DISTRIBUTION /AVAILABILITY  OF   REPORT 

Approved  for  public  release; 
distribution  is  unlimited 


a    PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


5    MONITORING  ORGANIZATION   REPORT  NUMBER(S) 


6a    NAME  OF  PERFORMING  ORGANIZATION 

Naval  Postgraduate  School 


6b    OFFICE  SYMBOL 
(If  applicable) 
Code    6  9 


7a    NAME  OF  MONITORING  ORGANIZATION 

Naval  Postgraduate  School 


6c    ADDRESS  (Ofy,  State,  and  ZIP  Code) 

Monterey,  California   93943-5000 


7b    ADDRESS  (Cry,  State,  and  ZIP  Code) 

Monterey,  California   93943-5000 


8a    NAME  OF  FUNDING  .  SPONSORING 
ORGANIZATION 


8b    OFFICE  SYMBOL 
(If  applicable) 


9    PROCUREMENT  INSTRUMENT  IDENTIFICATION   NUMBER 


8c.  ADDRESS  fOty,  State,  and  ZIP  Code) 


10    SOURCE  OF  FUNDING   NUMBERS 


PROGRAM 
ELEMENT  NO 


PROJECT 
NO 


TASK 
NO 


WORK    UNIT 
ACCESSION   NO 


11    TITLE  (Include  Security  Classification)  _  _     ,  „„    ,,    ^ 

AN    ADVANCED    STUDY    OF    NATURAL    CONVECTION    IMMERSION    COOLING    OF    A    3    x    3    ARRAY    OF 

SIMULATED  COMPONENTS  IN  AN  ENCLOSURE  FILLED  WITH  DIELECTRIC  LIQUID 


12,  PERSONAL  AUTHOR{S) 

Benedict,  Terry  J. 


13a.  TYPE  OF  REPORT 

Master's  Thesis 


13b    TIME  COVERED 
FROM  TO 


14    DATE  OF  REPORT   {Year,  Month,  Day) 

1988,    June 


15    PAGE  COUNT 
115 


16    SUPPLEMENTARY  NOTATION 

The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official 
policy  or  position  of  the  Department  of  Defense  or  the  U.S.  Government. 


17 


COSATI  CODES 


FIELD 


GROUP 


SUB-GROUP 


18    SUBJECT  TERMS  {Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Electronic  Cooling;  Protruding  Heat  Sources; 
Flow  Visualization;  Convective  Heat  Transfer; 
Immersion  Cooling;  Dielectric  Liquids 


19   ABSTRACT  {Continue  on  reverse  if  necessary  and  identify  by  block  number) 

An  experimental  study  has  been  conducted  to  examine  the  three-dimensional 
natural  convection  heat  transfer  from  an  array  of  simulated  electronic  com- 
ponents immersed  in  a  chamber  filled  with  Fluorinert  FC-75,  a  commercially 
available  dielectric  liquid.   The  top  and  bottom  walls  of  the  chamber  were 
maintained  at  uniform  temperature  while  all  other  surfaces  were  insulated. 
The  simulated  components  were  in  the  foinn  of  a  3  x  3  array  of  discrete  pro- 
truding aluminum  blocks,  each  with  geometrical  dimensions  of  a  20  Pin  Dual- 
inline-Package.   The  components  were  electrically  powered  resulting  in  a 
range  of  energy  dissipation  levels  from  0.1  to  3.1  watts.   Flow  visualization 
in  steady  state  was  accomplished  using  Magnesium  particles  illuminated  by  a 
Helium  Neon  laser  plane.   Component  surface  temperature  measurements  allowed 
determination  of  the  heat  transfer  characteristics.   Timewise  fluctuations 


20    DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
JP  UNCLASSIFIED/UNLIMITED       D  SAME  AS  RPT 


D  DTIC  USERS 


21     ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


22a    NAME  OF  RESPONSIBLE  INDIVIDUAL 

Prof.  Matthew  D.  Kelleher 


22b  TELEPHONE  f/nc/ude  Area  Code) 

(408)     646-2530 


22c    OFFICE^ SYM 

Code    69K, 


SOL 


DO  FORM  1473.  84  MAR 


83  APR  edition  may  be  useci  until  exhausted 
All  other  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF  ^HlS  ^AGE 

Q  U.S.  Government  Printing  0'*ice     1986—606-243 

UNCLASSIFIED 


UWUijAbbilr'iED 


SCCUniTY  CLAMiriCATION  OF  THIS  PAOE 


#19  -  ABSTRACT  -  (CONTINUED) 

Of  temperature  at  several  locations  were  measured  with 
increasing  power  levels. 


11  UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


Approved  for  public  release;  distribution  is  unlimited 

An  Advanced  Study  of 
Natural  Convection  Immersion  Cooling 
of  a  3  X  3  Array  of  Simulated  Components 
in  an  Enclosure  Filled  with  Dielectric  Liquid 

by 

Terry  J.  Benedict 

Lieutenant,  United  States  Navy 

B.S.,  United  States  Naval  Academy,  1982 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  ENGINEERING  SCIENCE 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
June  1988 


ABSTRACT 

An  experimental  study  has  been  conducted  to  examine  the 
three-dimensional  natural  convection  heat  transfer  from  an 
array  of  simulated  electronic  components  immersed  in  a 
chamber  filled  with  Fluorinert  FC-75,  a  commercially 
available  dielectric  liquid.  The  top  and  bottom  walls  of 
the  chamber  were  maintained  at  uniform  temperature  while  all 
other  surfaces  were  insulated.  The  simulated  components 
were  in  the  form  of  a  3  x  3  array  of  discrete  protruding 
aluminum  blocks,  each  with  geometrical  dimensions  of  a  20 
Pin  Dual-inline-Package.  The  components  were  electrically 
powered  resulting  in  a  range  of  energy  dissipation  levels 
from  0.1  to  3.1  watts.  Flow  visualization  in  steady  state 
was  accomplished  using  Magnesium  particles  illuminated  by  a 
Helium  Neon  laser  plane.  Component  surface  temperature 
measurements  allowed  determination  of  the  heat  transfer 
characteristics.  Timewise  fluctuations  of  temperature  at 
several  locations  were  measured  with  increasing  power 
levels. 
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I.   INTRODUCTION 

A.   STATEMENT  OF  PROBLEM 

As  microelectronic  components  continue  to  decrease  in 
size  and  increase  in  component  density,  the  problem  of  heat 
dissipation  becomes  increasingly  more  critical  in  the  design 
of  circuits  and  components.  A  five-order-of -magnitude 
increase  in  circuit  integration  in  the  last  25  years  has 
been  associated  with  successive  revolutions  in  components 
technology. 

With  the  recent  development  of  CMOS,  industry  has  seen 
an  almost  three-order-of-magnitude  reduction  in  size,  an 
order-of-magnitude  increase  in  the  characteristic  component 
dimensions  and  perhaps,  most  importantly,  a  precipitous  drop 
in  transition  switching  energy  from  10"^  J  in  1960  to  nearly 
10"-^-^  J  in  present  day  devices.  It  is  anticipated  that 
power  dissipation  of  100  W  in  a  10  million  component  chip 
could  be  encountered  by  1990.  This  would  result  in 
component  heat  fluxes  in  excess  of  100  W/cm^  [Bar-Cohen, 
1987,  pp.  159-175].  The  strong  requirement  to  keep 
components  at  low  temperature  levels,  100°C  to  125°C,  does 
exist,  since  for  every  2  0°C  decrease  in  components 
temperature,  the  components  failure  rates  are  cut  in  half 
[Oktay,  1986,  pp.  36-42]. 


B .   BACKGROUND 

Air  cooling  is  the  standard  design  in  today's  systems. 
Advances  in  air  cooling  technology  now  include  impingement 
cooling  and  combined  conduction  and  air  cooling  technology 
which  offers  even  higher  cooling  rates  than  natural  or 
forced  convection  air  cooling.  However,  air  as  a  heat 
transfer  medium  is  extremely  limited  not  only  in  its  ability 
to  transfer  heat  but  also  in  its  capacity  to  carry  the 
transferred  energy  with  a  reasonable  temperature  increase 
[Baker,  1972,  pp.  213-222]. 

Special  component  surfaces  have  been  investigated. 
Features  such  as  Turbulence  Promoting  Fins  (TPF)  have 
significantly  higher  cooling  enhancements  when  used  with 
special  gases  such  as  Helium  as  the  closed  cycle  coolant 
[Kishimoto  et  al.,  1983,  pp.  286-93].  By  doing  so,  a  300% 
greater  dissipation  rate  can  be  reached  than  in  an  open 
cycle  with  air  as  the  coolant. 

Forced  convection  in  liquids  offers  greater  heat 
transfer  coefficients  than  does  natural  convection  in 
liquids  or  air.  However,  forced  convection  introduces  two 
major  problems:  (i)  greater  amounts  of  expensive  coolant 
are  required  and  (ii)  pumps  and  external  plumbing  increase 
the  size  and  reduce  the  reliability  by  introducing  the 
possibility  of  vibration  related  failures  to  the  system 
[Baker,  1973,  pp.  163-173]. 


Boiling  can  produce  very  large  scale  heat  transfer 
coefficients,  but  it  also  introduces  major  problems:  (i) 
boiling  restricts  the  physical  design  of  the  system,  (ii)  if 
the  surface  is  very  smooth,  as  most  electronic  components 
are,  a  high  degree  of  superheat  may  be  required  so  that 
nucleation  can  occur.  The  occurrence  of  periodic  boiling 
due  to  the  high  local  heat  flux  may  result  in  severe  local 
turbulence  and  mechanical  damage  and  (iii)  boiling  may 
result  in  the  concentration  of  impurities  at  the  nucleation 
site  [Baker,  1973,  pp.  163-173]. 

C.   STUDIES  ON  NATURAL  CONVECTION  COOLING  OF  ELECTRONIC 
DEVICES 

1.   Numerical  Studies 

Of   the  different  heat  transfer  modes  mentioned 

above,  immersion  cooling  by  natural  convection  offers  not 

only  high  heat  transfer  potential,  but  also  a  practical 

means  of  controlling  circuit  temperatures.   However,  only 

limited  heat  transfer  data  are  available  for  immersion 

cooling  of  discrete,  protruding  heat  sources  [Liu  et  al., 

1987,  pp.  309-30].   Early  studies  of  natural  convection  heat 

transfer  were  motivated  by  applications  such  as  energy 

conservation  and  solar  energy  collection.    The  pertinent 

length  scales  were  several  order  of  magnitude  higher  then 

those   applicable   in   cooling   of   electronic   equipment. 

However,  recent  studies  have  focused  on  computational  and 


experimental   studies   of   natural   convection   cooling   in 
electronic  components. 

[Han  et  al.,  1984,  pp.  21-28]  have  used  numerical 
simulation  to  investigate  the  transient  laminar  natural 
convection  in  a  rectangular  enclosure  with  a  localized 
constant  heat  flux  source.  [Acharya  et  al.,  1983,  pp.  855- 
66]  have  investigated  the  natural  convection  phenomenon  in 
an  inclined  square  box  shaped  enclosure  which  has  energy 
sources  in  it.  [Jaluria,  1985,  pp.  883-92]  investigated  the 
interacting  natural  convection  flows  generated  by  flush, 
isolated  thermal  energy  sources,  such  as  electronic 
components,  located  on  a  vertical  adiabatic  surface. 
[Jaluria,  1982,  pp.  223-27]  also  analytically  studied  the 
natural  convection  flow  due  to  multiple  isolated  heated 
elements  located  on  a  vertical  surface.  [Torok,  1984,  pp. 
49-57]  used  the  finite  element  technique  to  predict 
component  and  board  temperatures  and  heat  fluxes. 

[Kuhn  et  al.,  1986,  pp.  55-62]  have  investigated  the 
transient  three-dimensional  natural  convective  flow  and  heat 
transfer  within  a  rectangular  box  following  a  sudden 
temperature  change  on  a  rectangular  portion  of  one  vertical 
wall.  They  found  that  when  two  elements  are  symmetrically 
placed  on  the  wall,  the  distance  between  them  has  only  a 
slight  effect  on  the  mean  heat  transfer  distribution. 
However,  it  does  have  a  significant  effect  on  the  local 
Nusselt   number.     [Shakerin   et   al.,   1986,   pp.   17-22] 


investigated  the  natural  convection  flow  adjacent  to  a 
heated  wall  with  single  and  repeated,  two-dimensional, 
rectangular  roughness  elements.  [Oosthuizen  et  al.,  1987, 
pp.  1-7]  studied  the  heat  transfer  rates  from  an  element  of 
square  cross-sectional  shape  held  at  a  uniform  temperature 
and  mounted  on  a  vertical  adiabatic  wall  of  a  tilted  square 
enclosure.  They  found  that  the  mean  heat  transfer  rate  from 
the  block  did  not  vary  continuously  with  the  angle  of 
inclination. 

[Liu  et  al.,  1987,  pp.  486-97]  have  conducted  a 
finite  difference  numerical  study  of  natural  convection  flow 
in  a  rectangular  enclosure  of  an  array  of  simulated 
components  that  are  mounted  on  a  vertical  wall.  The  results 
of  this  numerical  study  have  been  presented  for  enclosure 
widths  of  18  and  30  mm.   It  was  concluded  that: 

-  The  temperature  field  in  the  enclosure  is  characterized 
by  boundary  layer  regions  surrounding  individual 
components . 

-  There  are  only  small  interferences  among  the  components, 
especially  for  the  lower  rows. 

-  The  maximum  temperatures  on  the  components  surfaces  are 
located   on   the   upper   horizontal   surfaces   of   the 
components.   Lower  temperatures  have  been  found  on  the 
lower  horizontal  faces  of  the  components  in  the  bottom 
row. 

-  The  results  for  the  30  mm  width  are  essentially  the  same 
as  those  for  the  18  mm,  indicating  only  a  minor 
dependence  on  the  enclosure  width. 

[Liu  et  al.,  1987,  pp.  309-30]  have  also  observed  a  local 
oscillatory   surface   temperature   response   in   the   same 


enclosure.    Their  study  found  that  the  components  surface 
temperatures  have  a  tendency  to  fluctuate  within  a  range  of 
approximately  3°C  with  a  period  of  4  seconds. 
2 .   Experimental  Studies 

Numerous  experimental  studies  have  been  conducted 
concurrently  with  the  numerical  studies  described  above. 
Two  of  the  earliest  studies  were  conducted  by  [Baker,  1972, 
pp.  213-22]  and  [Baker,  1973,  pp.  163-73].  The  purpose  of 
these  investigations  was  to  examine  liquid  immersion 
techniques  for  cooling  miniature  heat  sources,  specifically 
effects  of  component  size  on  the  convective  heat  transfer. 
[Goel  et  al.,  1986,  pp.  1341-46]  experimentally  investigated 
the  two-dimensional  natural  convection  flow  from  a  long, 
finite  size  heat  source  located  on  a  vertical  or  inclined 
plate  in  air.  The  principal  items  of  interest  were  the 
nonboundary  layer  effects  that  arise  in  the  flow, 
particularly  in  the  vicinity  of  the  source.  Flow  separation 
was  observed  at  large  inclinations. 

[Milanez  et  al.,  1986,  pp.  1347-52]  experimentally 
studied  the  interaction  of  two  horizontal,  wide  strips  on  a 
vertical  adiabatic  wall  in  both  air  and  water.  They  found 
the  results  of  their  experiment  in  good  agreement  with  past 
numerical  studies.  [Ortega  et  al.,  1986,  pp.  5-15] 
presented  the  results  for  a  series  of  experiments  in  air 
designed  to  measure  the  natural  convective  heat  transfer 
rates  and  the  thermal  characteristics  of  a  rectangular  array 


of  cubical  roughness  elements  mounted  on  an  insulated  plane 
wall,  with  and  without  a  shrouding  wall.  [Parks  et  al., 
1987,  pp.  91-109]  simulated  microelectronic  circuits  with 
thin  foil  heaters  for  both  flush  and  protruding  components. 
They  found  that  the  protruding  heaters  have  approximately 
15%  higher  coefficients  and  that  the  upper  heaters  have 
varying  coefficients  depending  upon  the  array  configuration. 
A  number  of  studies  have  been  conducted  at  the  Naval 
Postgraduate  school.  [Knock,  1983]  studied  the  effect  of 
the  location  of  a  single  protruding  heater  in  an  enclosure 
filled  with  water.  [Pamuk,  1987]  investigated  the  heat 
transfer  characteristics  of  an  array  of  nine  protruding 
components  immersed  in  a  dielectric  fluid  in  an  uninsulated 
chamber.  [Hazard,  1987]  studied  the  natural  convection 
liquid  cooling  of  a  vertically  arranged  components  array 
with  and  without  a  shrouding  wall. 

D.   OBJECTIVES 

This  experimental  work  was  motivated  by  the  numerical 
study  of  [Liu  et  al.,  1987,  pp.  486-97]  and  the  recent 
experimental  studies  of  [Pamuk,  1987].  The  objectives  for 
this  thesis  are: 

-  To  design  and  build  a  circuit  card  of  nine  simulated 
components  and  determine  the  rate  of  heat  transfer  for 
various  power  inputs. 

-  To  determine  the  natural  convection  flow  patterns 
developed  in  the  chamber  by  means  of  flow  visualization. 

-  To  correlate  the  data  in  terms  of  dimensionless  heat 
transfer  parameters  such  as  Rayleigh  vs.  Nusselt  number. 


II.   EXPERIMENTAL  APPARATUS 

A.   GENERAL  CONSIDERATIONS 

This  experimental  study  is  a  continuation  of  past 
experiments  conducted  at  the  Naval  Postgraduate  School 
[Knock,  1983;  Pamuk,  1987;  Hazard,  1987].  The  purpose  of 
this  work  was  to  examine  natural  convection  heat  transfer 
and  fluid  flow  from  a  3  x  3  element  array  of  simulated 
electronic  components.  A  major  change  from  the  previous 
experiments  was  the  design  of  the  simulated  circuit  card. 
In  [Pamuk,  1987]  an  arrangement  of  vertical  components  was 
used.  The  present  study  utilized  a  horizontal  arrangement. 
A  decision  to  increase  the  card  thickness  was  made  to  reduce 
the  contamination  of  the  dielectric  fluid  by  the  loose 
insulating  material  previously  used  and  also  to  reduce  the 
conduction  loss  through  the  card. 

The  experimental  apparatus  consists  of  a  chamber  filled 
with  a  dielectric  fluid.  Top  view  of  the  chamber  is 
provided  in  Figure  2.1.  The  top  and  bottom  surfaces  of  the 
chamber  are  maintained  at  a  prescribed  temperature  by  means 
of  two  heat  exchangers  and  two  separate  chilled  water 
circulation  baths.  One  vertical  wall  of  the  chamber 
simulates  a  circuit  card.  It  contains  nine  simulated 
components  in  a  3  x  3  array,  while  all  other  vertical  walls 
are  insulated.   The  geometrical  arrangement  of  the  simulated 


1 

^ ^E-iJ 

S        1   21 

-1           r 

r- 

'^" 

UJ 

o 

< 

CO 

CO 

^ 

< 

Q. 

O 

z 

cr 

$ 

E 
E 

/ 

-^ 

T 

-. 

^ 

V 

/ 

1^. 

E 
E 

CNJ 

1 
O 

o 

E 

cr 

E 

< 

<o 

o 

■ 

" 

-I             u 

i_ 

f                  1 

^                                E 
to 

•w 

u 

0) 


Id 
xi 
o 

o 

(U 
•H 
> 

o 

Eh 


(U 

p 

en 

•H 

[X4 


components  on  the  card  is  seen  in  Figure  2.2.  The 
components  were  chosen  to  simulate  an  array  of  2  0  pin  Dual- 
inline-packages  (DIPS) .  Their  dimensions  are  identical  to 
those  used  by  [Liu  et  al.,  1987,  pp.  486-97;  Liu  et  al., 
1987,  pp.  309-30;  Pamuk,  1987]  in  order  to  permit  future 
comparisons  of  experiments  and  computations. 

The  heating  elements  were  made  of  aluminum  and  heated  by 
means  of  small  foil  type  resistance  heaters  attached  to  the 
bottom  surfaces.  The  temperatures  of  the  faces  of  the 
elements  were  measured  by  means  of  embedded  thermocouples. 
Data  were  recorded  and  processed  using  an  automatic  data 
acquisition  system.  The  major  components  of  the  system  are 
described  in  the  following  section. 

B .   COMPONENTS 

1.  Heating  Element 

Each  heating  element  was  a  block  of  aluminum  6mm  by 
8mm  by  24mm  (Figure  2.3).  The  temperatures  at  the  faces  of 
the  blocks  were  measured  by  means  of  thermocouples  placed 
within  grooves  (0.5mm  deep)  on  each  surface.  The  front 
surface  thermocouple  was  passed  through  the  block  by  means 
of  a  1.02mm  hole  drilled  through  the  block.  All  grooves 
were  filled  with  high  thermal  conductivity  epoxy  and 
smoothed  prior  to  attachment  to  the  simulated  circuit  card. 

2 .  Thermocouples 

All  thermocouples  were  10-mil  (0.254mm)  copper- 
constantan.   These  small  wires  were  chosen  to  minimize  the 
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Figure  2.3   Heating  Element 


12 


material  removed  from  the  block.  The  thermocouples  were 
manufactured  by  using  an  arc  welder  designed  for  such 
purposes.  The  setting  for  both  current  time  and  gas 
pressure  (argon)  was  kept  at  'medium. •  Due  to  the  extremely 
small  diameter  of  the  wire  the  gas  pressure  was  carefully 
regulated  to  less  than  5  psig.  The  thermocouples  were  also 
checked  for  continuity  due  to  the  fact  that  the  arc  welder 
tended  to  embrittle  the  copper  wire  at  the  bead. 

3 .  Heaters 

A  nearly  constant  heat  flux  condition  was  maintained 
at  the  base  of  the  block  by  use  of  thin  (0.18mm)  foil  type 
heaters  (Figure  2.4).  Each  of  the  heaters  was  measured  to 
have  a  resistance  between  10.5  and  10.6  ohms.  These  heaters 
were  designed  to  align  with  the  blocks.  Holes  were  provided 
within  the  heaters  to  pass  the  various  thermocouple  wires. 

The  power  leads  of  the  heaters  were  gold  plated  to 
ensure  a  low  electrical  contact  resistance.  Wires  from  the 
power  supply  were  attached  to  the  leads  by  means  of 
soldering.  After  ensuring  proper  alignment,  the  heaters 
were  bonded  to  the  blocks  by  means  of  a  high  thermal 
conductivity  epoxy. 

4.  Simulated  Circuit  Card 

Three  grooves  of  0.5mm  depth  were  made  on  the  front 
of  the  card  in  order  to  allow  the  heaters  to  be  mounted 
below  the  surface  and  thus  prevent  the  dielectric  fluid  from 
seeping  under  the  blocks  and  introducing  a  variable,  unknown 
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resistance.   The  heating  elements  were  bonded  to  the  circuit 
card  with  a  low  thermal  conductivity  epoxy. 

5.  Test  Chamber 

The  test  chamber  was  constructed  of  plexiglass  (13mm 
thick)  and  had  inside  dimensions  of  120mm  in  height,  144mm 
in  length,  and  3  0mm  in  width.  The  wires  were  fed  through 
the  grooves  in  the  back  of  the  circuit  card  to  a  "Tygon" 
tube  on  the  back  wall.  The  chamber  also  contained  two  small 
pressure  equalization  outlets,  located  at  the  top  and  the 
bottom  of  the  front  vertical  surface. 

6.  Heat  Exchangers 

The  two  aluminum  heat  exchangers  measure  3  8mm  by 
65min  by  274mm.  The  details  of  the  heat  exchangers  are 
available  in  [Knock,  1983].  Aluminum  plates  3mm  thick  were 
used  as  top  and  bottom  surfaces.  These  plates  were  grooved 
2.5mm  to  provide  slots  for  the  thermocouples.  This  left  a 
0.5mm  wall  thickness  of  aluminum  between  the  thermocouple 
bead  and  the  surface  in  contact  with  the  dielectric  fluid. 
The  resulting  conductive  thermal  resistance  was  found  to  be 
negligible.  The  heat  exchangers  were  each  attached  to 
individual  coolant  baths  by  means  of  "Tygon"  tubing.  Four 
10  mil  (0.2  54mm)  thermocouples,  two  on  the  top  and  two  on 
the  bottom,  were  used  to  read  the  heat  exchanger 
temperatures . 
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7.  Flow  Visualization 

Visualization  was  accomplished  with  magnesium 
powder,  325  mesh.  Magnesium  has  a  density  of  1.92  q/crn^ 
compared  to  1.76  g/cm-^  for  FC-75,  resulting  in  a  large 
suspension  time  of  the  particles.  The  beams  from  the  light 
sources  were  spread  into  plane  sheets  which  illuminated  the 
particles  suspended  in  the  FC-75.  This  technique  allowed 
for  the  visualization  of  a  single  plane  of  the  flow  field. 
Other  planes  were  examined  by  the  realignment  of  the 
apparatus.  A  four  milli-watt  helium  neon  laser  and  a 
cylindrical  lens  were  used  for  the  time  lapse  camera  photos 
(Figure  2.5) . 

8 .  Assembly 

After  sliding  the  card  into  the  slot  in  the  chamber, 
both  the  thermocouple  and  power  wires  were  run  out  through 
the  1.27cm  diameter  extension  pipe  on  the  back  of  the 
chamber  to  the  connected  "Tygon"  tubing  of  approximately  1.0 
meter  in  length. 

Rubber  o-rings  were  placed  in  grooves  located  on  the 
top  and  bottom  surfaces  of  the  chamber  and  the  thin  aluminum 
plates  were  set  in  place.  The  heat  exchangers  were  then  set 
on  the  plates  and  clamped  by  means  of  specially  made 
aluminum  clamps  to  ensure  a  proper  seal.  The  inlet  and 
outlet  tygon  tubes  were  attached  to  the  constant  temperature 
baths   for  each  heat   exchanger.   Finally,   the   pressure 
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equalization   tubes   were   attached   to   the   chamber.     a 
schematic  of  the  entire  assembly  is  shown  in  Figure  2.6. 
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III.   EXPERIMENTAL  PROCEDURE 

A.   APPARATUS  PREPARATION 

After  assembling  the  apparatus,  the  chamber  was  filled 
with  the  dielectric  fluid  through  the  top  vent  tube  by  means 
of  a  small  funnel  and  a  hypodermic  syringe.  To  damp 
disturbances,  the  entire  apparatus  was  placed  on  a  separate 
platform.  This  platform  also  served  as  a  means  of  leveling 
the  chamber  by  use  of  adjustable  screw  legs.  Since  the 
entire  apparatus  was  sealed  by  means  of  four  clamps  it  could 
be  tilted  in  order  to  remove  the  last  remaining  air  bubbles 
inside  the  chamber. 

The  use  of  a  single  constant  temperature  bath  providing 
cooling  water  in  series  to  the  top  and  bottom  heat 
exchangers  proved  inefficient  to  maintain  both  the  top  and 
the  bottom  heat  exchangers  at  the  same  constant  temperature. 
A  second  bath  was  added  so  that  each  heat  exchanger  was 
supplied  individually  and  could  be  adjusted  separately.  The 
ten  terminals  on  the  power  supply  panel  were  used  to  supply 
power  to  the  heaters  on  the  simulated  components.  The  first 
terminal  was  the  main  input,  the  next  nine  were  the 
terminals  to  the  individual  heaters  located  on  each 
component.  These  heaters  were  connected  in  series  with  the 
2  ohm  precision  resistors.   This  arrangement  was  repeated 
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for  each  component.   These  nine  units  were  then  connected  in 
parallel  to  the  main  input. 

B.   DATA  ACQUISITION 

Prior  to  energizing  the  power  supply  or  the  constant 
temperature  baths  all  the  thermocouples  were  scanned  to 
measure  the  starting  steady  state  temperatures.  The 
constant  temperature  baths  were  set  at  10°C  by  means  of 
adjustment  dials  located  on  the  baths.  Throughout  each 
experiment  the  bath  supplying  the  top  heat  exchanger 
required  constant  monitoring.  As  the  power  levels 
increased,  the  bath  temperature  was  reduced  accordingly  to 
maintain  the  top  heat  exchanger  thermocouple  readings  at  a 
constant  10°C. 

The  desired  voltage  was  obtained  by  setting  the  power 
supply  to  the  desired  level.  The  reading  on  the  power 
supply  was  verified  by  means  of  reading  channel  number  61, 
the  main  power  input  channel,  on  the  data  acquisition  panel. 

Approximately  four  hours  were  allotted  for  the  system  to 
achieve  steady  state.  For  this  experiment,  the  steady  state 
condition  was  defined  as  a  majority  of  the  thermocouple 
readings  changing  less  than  0.2°C  during  three  separate 
readings  taken  at  ten  minute  intervals.  The  bath 
temperatures  were  checked  to  ensure  they  were  at  10°C.  If 
they  were  in  agreement  then  the  data  was  recorded  by  means 
of  the  program  "ACQUIRE"  (Appendix  C) .  This  program  scanned 
all  76  channels  and  converted  the  emf  signals  from  the 

21 


thermocouples  to  temperatures  in  degrees  Celsius.  The  power 
to  the  heater  was  determined  by  utilizing  the  precision  2 
ohm  resistor.  This  resistor  was  connected  in  series  with 
the  heater  and  allowed  the  measurement  of  the  current 
through  this  configuration.  This  current  was  then 
multiplied  by  the  measured  voltage  of  the  heater  to  produce 
the  power  input  to  the  component. 

A  second  program  "Fastscan"  (Appendix  C)  was  developed 
to  scan  three  individual  channels.  Each  of  the  three 
channels,  numbers  1,  13,  32,  was  scanned  individually  at  a 
rate  of  3  scans  per  second  for  a  period  of  200  seconds.  The 
main  intention  in  doing  so  was  to  isolate  any  periodic 
behavior  caused  by  the  flow  past  these  three  thermocouples. 
A  third  program  "Plot"  (Appendix  C)  was  developed  to  plot 
these  results  on  a  HP  model  7470A  plotter. 

Flow  visualization  was  conducted  by  means  of  325  mesh 
Magnesium  particles  and  a  four  milliwatt  laser  focused  with 
a  cylindrical  lens.  The  magnesium  particles  were  injected 
into  the  chamber  by  way  of  the  bottom  vent  tube.  A  motor 
driven  Nikon  F3  camera  with  a  Micro-Nikkor  55mm  lens  and  a 
Nikon  MT-2  Intervalometer  was  used  for  the  long-time 
photography.  The  time  exposure  photographs  of  length  8  or 
2  0  seconds  were  taken  at  the  various  power  levels  where  heat 
transfer  data  had  been  obtained.  The  photographs  were  taken 
of  the  three  columns  of  components  by  placing  the  camera  at 
a  ninety  degree  angle  to  the  plane  created  by  the  laser. 
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This  configuration  was  then  reversed  to  obtain  photographs 
of  three  separate  planes  in  the  chamber. 

C.   DATA  ANALYSIS 

1.   DeterTnination  of  the  Nusselt  Number 

The  Nusselt  Number  is  defined  as  the  following: 


Nu  =  ^  (3.1) 

^f 


where  L  is  a  characteristic  length,  taken  here  as  the  length 
of  the  component  in  the  vertical  direction,  kf  is  the 
thermal  conductivity  of  the  fluid,  and  h  is  the  convective 
heat  transfer  coefficient.  All  fluid  properties  were  taken 
to  be  a  function  of  the  average  film  temperature. 

To  evaluate  the  heat  loss,  one  thermocouple  was 
placed  directly  behind  each  of  the  heated  blocks.  A  series 
of  nine  thermocouples  were  also  attached  to  the  back  of  the 
card.  These  temperatures  were  used  to  calculate  the 
conduction  loss  through  the  back  of  the  card. 

One  dimensional  conduction  resistance,  Rcon'  ^°^  ^^® 
simulated  board  was  determined  from: 


^con  ~  kA  (3.2) 


where  L  is  the  length  along  the  conduction  path  i.e.,  the 
thickness  of  the  board,  k  is  the  thermal  conductivity  of  the 
material   (Appendix  A)   and  A  is  the  area  normal  to  the 
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conduction  path.   The  area  normal  to  the  path  of  conduction 

was  calculated  as  the  area  directly  behind  the  components 

back,  192inin^.   By  doing  so  spreading  losses  were  neglected. 
The  conduction  loss  is  then: 


T   -  T, 
s    b 

^con  =  ^ (3.3) 

con 


Power  input  to  each  component  is  calculated  from: 


V.   -  V,  V, 
>.   -    in     h   h 
!in  -  R (3.4 


where  Vj^^  i^  the  input  voltage  applied  to  both  the  2  ohm 
resistor  and  the  heater  which  is  in  series  with  it.  Vj^  is 
the  voltage  drop  across  the  heater.  The  difference  between 
V-Ln  ^^d  ^h  gives  the  voltage  drop  across  the  2  ohm  resistor. 
The  net  energy  added  to  the  fluid  in  the  chamber  is 
the  difference  between  Qj^^  ^^^  Qcon  which  by  Newton's  Law  of 
cooling  is: 

Qconv  ~  Qnet  ~  ^^C^avg  ~  "^c) 

=  Qin  -  Qcon  (3.5) 


where, 


A  =  The  total  wetted  surface  area  of  the 
components , 

'avg  ~  Area  weighted  average  surface  temperature. 
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calculated  by, 


avg 


E  A.T. 
1  1 

L    A. 


(3.6) 


where , 


rp.   = 


A^  = 


Temperature  of  the  i^^  surface  of  the 
components, 

Area  of  the  i^^  surface  of  the  components 
The  summation  of  the  index  i  varies  from 
1  to  5, 


and. 


Thus, 


Tq  =   Average  heat  exchanger  temperature. 


h  = 


'conv 


A.  .  (T     -  T^) 
tot   avg     c 


(3.6) 


The  above  procedure  was  repeated  for  each  of  the 
nine  components  at  each  power  level.  The  Nusselt  Number  was 
then  calculated  from  Equation  3.1. 

2 .   Determination  of  the  Rayleiqh  Numbers 

The  definition  of  the  temperature  based  Rayleigh 
number  is: 


Rat  =  Grt  *  Pr 


(3.7) 
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where  the  Grashof  Number  (Gr^)  is  defined  as, 

g3L^(T     -  T  ) 
Grt  =  ^^^2 2_  (3.8) 


and  the  Prandtl  Number  (Pr)  is  defined  as, 


Pr  =  ^  (3.9) 


where. 


a   = 


(pcp) 


All   the  properties   in  the   above   equations   are 
calculated  at  the  average  film  temperature: 


T     +  T 
Tf  =   ^^g  ,   ^  (3.10) 


Therefore: 


g3L^ (T     -  T  ) 
Rat  =  ^^ —  (3.11) 


3 .   Determination  of  the  Flux  Based  Rayleiqh  Number 

The  definition  of  the  flux  based  Rayleigh  number  is: 

Raf  =  Grf  *  Pr  (3.12) 
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all  film  properties  are  calculated  at  the  average  film 
temperature  as  before  and: 

Grf  =  2 '   (3.13) 

f    tot 

Functional  relationships  for  the  properties  at  the 
film  temperature  are  given  in  Appendix  A,  together  with 
sample  calculations.  Uncertainty  calculations  are  given  in 
Appendix  B. 
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IV.   RESULTS 

A.   FLOW  VISUALIZATION 

Flow  visualization  was  accomplished  using  Magnesium 
particles  as  a  flow  marker  and  a  Helium  Neon  laser  to 
illuminate  the  particles.  The  results  provided  insights 
into  the  flow  patterns  existing  at  various  power  dissipation 
levels.  As  expected,  the  particle  traces  suggest  highly 
three-dimensional  flow  patterns  at  all  power  levels. 
Particles  can  be  seen  entering  and  leaving  the  plane  of 
light  in  almost  all  the  long  time  exposure  photographs.  A 
multitude  of  flow  patterns  was  observed  for  the  wide  range 
of  input  power  levels  in  the  of  0.1  to  3.1  watts. 
Representative  flow  visualizations  are  presented  in  Figures 
4.1  to  4.14. 

1.   Flow  Patterns  at  Component  Midplanes 

The  following  flow  patterns  are  in  a  vertical  plane 
through  the  center  of  components  4,  5,  6  (Fig.  2.2).  The 
laser  was  placed  in  front  of  the  chamber  and  the  camera 
located  to  the  side  (Fig.  2.5).  The  camera  exposure  time 
was  set  at  10  seconds.  The  three  dark  regions  located  on 
the  left  side  of  the  photographs  are  the  three  rows  of 
components  as  seen  in  side  view.  The  top  and  bottom  of  the 
chamber  and  the  horizontal  walls  can  also  be  seen  in  the 
photos. 
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Figure  4.1  is  the  flow  pattern  for  a  power  input  of 
0.1  watts  and  indicates  a  pattern  that  contains  two 
individual  eddies  located  between  the  rows  of  components 
along  with  an  overall  circular  pattern  that  encompasses  the 
entire  chamber.  A  boundary  layer  like  flow  is  seen  in  the 
vicinity  of  the  heated  components.  A  dark  region  located  at 
the  bottom  of  the  photograph  is  evident  in  all  photos  taken 
throughout  the  experiment.  This  region  is  the  stagnant, 
dense  cold  fluid  layer  in  the  vicinity  of  the  bottom  heat 
exchanger.  This  region  proved  to  be  impermeable  at  all 
power  levels. 

Figure  4.2  for  power  input  of  0.3  watts  indicates 
the  same  overall  pattern  with  the  exception  of  a 
constriction  of  the  overall  circular  pattern  at  the 
approximate  height  of  the  second  row  of  components.  The 
slight  increase  in  the  upper  eddy  size  is  also  evident. 

Figure  4.3  for  a  power  input  of  1.5  watts  indicates 
an  increase  in  the  thickness  of  the  boundary  layer  emanating 
from  the  second  row  of  components  along  with  increased 
entrainment  of  the  particles  located  below  the  second  row  of 
components.  Also  evident  is  the  disappearance  of  the  two 
distinct  eddies  seen  earlier  and  strong  evidence  of  three 
dimensionality  in  the  flow  as  seen  by  the  crossing  of  the 
particle  paths. 
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Figure  4.1   Power  §  0.1  W 
Side  View 


Figure  4.2   Power  @  0.3  W 
Side  View 
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Figure  4 . 3 


Power  @  1 . 5  W 
Side  View 
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Figure  4.4,  for  a  power  input  of  2.7  watts,  shows  a 
buoyant  fluid  layer  adjacent  to  the  components  and  random 
motion  throughout  the  remaining  chamber.  Also  seen  are 
entrainment  of  particles  at  the  approximate  height  of  the 
second  row  of  components,  and  evidence  of  an  area, 
immediately  above  the  black  region  at  the  photo's  bottom, 
where  particles  are  basically  suspended. 

Figure  4.5,  at  3.1  watts,  demonstrated  a  downward 
flow  pattern  in  the  vicinity  of  the  bottom  row.  This 
downward  flow  was  forced  past  the  face  of  the  bottom  chip  in 
this  column.  The  flow  also  demonstrated  a  perpendicular 
entrainment  of  particles  between  the  first  and  second  row  of 
components,  and  a  large  increase  in  the  boundary  layer  at 
the  third  row  of  components. 

Figures  4 . 6  and  4 . 7  show  flow  patterns  in  the 
vertical  plane  at  the  mid-span  of  components  7,  8  and  9. 
These  two  photographs  were  taken  approximately  5  seconds 
apart  at  a  power  setting  of  3 . 1  watts.  They  demonstrate  the 
strongly  time-dependent  behavior  of  the  flow.  In  Figure 
4.6,  entrainment  into  the  upward  moving  buoyant  layer  is 
seen  to  be  present  almost  as  far  down  as  the  bottom  of  the 
lowest  component.  Far  away  from  the  components  the  core 
fluid  exhibits  no  preferential  direction  of  motion.  In 
Figure  4.7  the  flow  pattern  has  changed.  Near  the  lowest 
row  of  components  the  flow  now  appears  almost  quiescent. 
Also,  the  buoyant  region  appears  thinner  than  in  Figure  4.6. 
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Figure  4.4   Power  @  2.7  W 
Side  View 


Figure  4 . 5   Power  Q  3 . 1  W 
Side  View 
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Figure  4.6 


Power  @  3 . 1  W 
Side  View, 
Chips  7,  8,  9 


Figure  4 . 7 


Power  @  3 . 1  W 
Side  View, 
Chips  7,  8,  9 
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2 .   Flow  Patterns  in  Other  Planes 

Figure  4.8  is  a  schematic  of  the  three  planes  that 
were  investigated  when  the  positions  of  the  camera  and  the 
laser  were  reversed.  Figures  4.9,  4.10  and  4.11  are  at  a 
power  setting  of  0.1  watts  and  demonstrate  the  flow  patterns 
in  the  back,  middle  and  front  planes,  respectively.  In 
Figure  4.9,  the  back  plane,  the  flow  appears  to  be  somewhat 
like  a  large  circular  pattern  that  is  bounded  by  the 
physical  dimensions  of  the  chamber  and  exhibits  no 
noticeable  flow  directly  up  the  center  of  the  chamber. 
Figure  4.10,  the  middle  plane,  indicates  a  strong  upward 
flow  in  the  center  of  the  chamber  and  two  areas  of  eddy 
formation  located  approximately  mid-height  in  the  chamber. 
Finally,  Figure  4.11,  the  front  plane,  demonstrated 
approximately  the  same  behavior  as  the  middle  plane;  a 
general  area  of  motion  in  the  center  of  the  chamber  along 
with  eddy  formation  areas  on  either  side. 


^ 
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Figure  4.8   Top  View  of  Planes 
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Figure  4.9   Back  Plane  @  0.1  Watts 


36 


Figure  4.10   Middle  Plane  @  0.1  Watts 
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Figure  4.11   Front  Plane  @  0.1  Watts 
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Figures  4.12,  4.13  and  4.14  are  the  same  three 
planes  at  a  power  setting  of  3.1  watts.  Figure  4.12,  the 
back  plane,  indicates  regions  of  upward  flow  near  the 
components.  These  patterns  characterize  the  boundary  layer 
nature  of  transport  in  the  vicinity  of  the  elements. 
Relatively  short  particle  traces  in  regions  away  from  the 
components  appear  to  result  due  to  the  strongly  three- 
dimensional  transport.  The  dark  area  in  the  very  bottom  of 
the  photograph  is  again  pointed  out  to  demonstrate  the  cold 
isothermal  layer  developed  by  the  bottom  heat  exchanger. 

Figures  4.13  and  4.14  are  views  of  the  center  and 
front  planes  respectively.  The  flow  in  both  planes  appears 
strongly  time  dependent.  Figure  4.14  shows  relatively 
longer  particle  traces  than  Figure  4.13,  indicating  stronger 
three  dimensional  effects  near  the  components. 

These  13  photographs  illustrate  the  complex  three- 
dimensional,  time  dependent  flow  patterns  found  in  this 
study.  Obtaining  a  detailed  understanding  of  these  effects 
is  essential  in  achieving  predictive  capability  for  the  heat 
transfer  in  these  configurations. 

B.   REDUCTION  OF  HEAT  TRANSFER  DATA 

The  temperature  and  reduced  data  tables  corresponding  to 
the  eight  power  levels  are  as  follows:  (see  Figure  2.2  for 
identification  of  the  component  numbers) . 
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Figure  4.12   Back  Plane  @  3.1  Watts 
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Figure  4.13   Middle  Plane  @  3.1  Watts 
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Figure  4 . 14   Front  Plane  @  3 . 1  Watts 
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TABLE  1 
TEMPERATURE  DATA  FOR  INPUT  POWER  Q^j^  =  0.3  W 

■r^ZZZ    ^E:uLI5  riftE  ITCRED  IN  FILE:   D809Bi;:0 


r^r^PIENT  ''E^P  UAS:   ^3  c 
VCL'.'^E'E^?  ^E^DIN'.  -AS:   J. 0382  «hl 
BAiH  TTMp  ^^c..    ,,)  ,- 

THi:  RUN  knS  RECORDED  UITH  TLOU  VIZ 

ALL  TEMPERAFURES  ARE  IM  DEGREES  CELCIUS 

CENTER     TOP         RIGHT      LEFT 

CHIPNOI:  IS.S'JE^OO   l6.33E*nO   1S.I5E*00   16.1bE*00 
POWER  (WATTS):   23.33E-02 

CHIP  N02:  I6.58E»00   t6.31E*00   I5.34E+00   16.19E*00 
POWER  (WATTS):   29.3^E-02 

CHIP  N03:  16.74E*00   16.:2E*00   I6.^6E+00   I6.33E*00 
POWER  (WATTS):   23.83E-02 

CHIP  N04:  1S.E2E*-00   IS.33E*00   16.03E*00   15.82E*00 

POWER  (WATTS):   23.55E-02 

CHIP  NOS:  1S.52E*00   tG.;SE*00   16.!2E*00   !G.05E*00 
POWER  (WATTS):   23.b7E-02 

CHIP  N06:  I6.5IE*00   I6.50E>00   1S.33E^00   18.33E-00 

POWER  (WATTS):   29.e^£-02 

CHIP  N07:  1S.S4E*00   16.42E*00   I6,I2E*00   16.12E*00 
POWER  (WATTS):   30.01E-02 

CHIP  N08:  18.71E*00   1G.4:E*00   t6.40E*00   1G.40E*00 
POWER  (WATTS):   30.12E-02 

CHIP  N03:  lS.bSE*00   IG.3SE*00   t6.23E*00   16.3GE*00 
POWER  (WATTS):   29.37E-02 


BOTTOM 

BAC^^ 

16. I 1E*00 

16 

.  '36£*00 

16.30E*00 

17 

.53E*00 

16.63E+00 

17 

.40E*00 

IG.02E^OO 

16 

.33E*00 

16.30E*00 

17 

.  i  ;£*oo 

1G.58£*00 

18 

13E>00 

16.32E*00 

17 

3!E*00 

ie.28E*00 

18 

02E*00 

16.41E*00 

17 

63E*00 

HEAT  EVCHANGERS  TEMPERATURES; 

bottom:   • 
TCP: 


RIGHT      LEFT 
99.12E-01   99.39E-01 
10.33E*00   98.a3E-01 


BACK  PLANE  TEMPERATURES  ARE 


T(55): 

14 

67E*00 

T(56): 

15 

.^1E*00 

T(57): 

15 

27E+00 

T(72): 

15 

57E+00 

T(73): 

15 

55E»00 

T(74): 

16 

06E*00 

T(75): 

15 

56E*00 

T(76): 

15 

56E*00 

T(77): 

1  c 

aSE  +  ilO 
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TABLE    2 
TEMPERATURE    DATA    FOR    INPUT    POWER   Q^^    =    0.7    W 

THESE   RESULTS   ARE   3T0RED    IN   FILE:      D30982:00 

AMBIENT    TEhP   UAS :         23   C 
VCLT.iETER   RErtOKiG   WAS:      3. 2198    *bl 
BATH    TE^^P    UAS:  '0    C 

THi:    RUN    WAS    RECORDED    WITH    FLOW    VIZ 

ALL    TEMPERATURES   ARE    IN   DEGREES    CELCIUS 

CENTER  TOP    '  RIGHT  LEFT  BOTTOM  BACK 

CHIP   N01:    23.35E*00      22.53E*00      22.35E*00      22.39E>00      22.24E*00  2A.06E*00 

POWER    (WhTTS):      71.08E-02 

CHIP   N02:    23.36E*00      22.37E*00     22.43E*00      22.09E*00      22.39E*00  25.05E+00 

POWER    (WATTS):      71.11E-02 

CHIP   N03:    23.I6E*00      22.3GE»00      22.55E*00      22.ME»00     22. 73E*00  24.5aE  +  00 

POWER    (WATTS):      70.87E-02 

CHIP  N04:    23.23E*00     22.42E*00     22.31E*00      21. 69E-00      22.IGE*00  24.16E+00 

POWER    (WATTS) :      70.18E-02 

CHIP   NOS:    23.nE*00      21.96E*00      21.94E  +  00      21.8SE*00      22.37E  +  00  24.igE*00 

POWER    (WATTS):      70.47E-02 

CHIP  N06:   22. S9E*00     22.77E*00     22.18E*00      22.21E*00     22.SOE*00  26.12E+00 

POWER    (WATTS):      70.38E-02 

CHIP   N07:    23.49E*00     22.7SE*00     22.I7E*00      22.15E*00     22.39E*00  24.55E*00 

POWER    (WATTS):      71.26E-02 

CHIP  N08:  23.32E*00  22.G3E>00  22.47E+00  22.55E*00  22.47E>00    26.17E*00 
POWER  (WATTS):   71.55E-02 

CHIP  N09:  23.07E-00  22.5'iE*00  22.05E>00  22.38E*00  22.20E*00    25.4IE*00 
POWER  (WATTS):   70.35E-02 

HEAT  EXCHANGERS  TEMPERATURES:     RIGHT      LEFT 

BOTTOM:  98.18E-01   98.A6E-01 

TOP:     _  98.21E-0I   93.S8E-0t 

BACK  PLANE  TEMPERATURES  ARE  : 

T(55):  17'.67E*00 

T(55):  18.76E*00 

T(57):  t9.05E*00 

T(72):  19.08E*00 

T(73):  18.90E*00 

T(74):  19.nE*00 

T(75):  t9.02E*00 

T(76):  18.79E*00  " 

T(77):  18.7tE*00 
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TABLE    3 
TEMPERATURE    DATA    FOR    INPUT    POWER    Q^^    =    1.1    W 


ThESE    RESULTS   AR£   STORED    IN   FILE:      D8099074S 

AMBIENT    TE«P   wfiS:        23   C 
vOlI.IETE^   '^"E'^DING   UAS:      ^.1107    *hl 
BATH    TEr-P    UAS:  IQ    r 

THIS   RUN   WAS   RECORDED   WITH   FLOW   VIZ 

ALL    TEMPERATURES  ARE    IN   DEGREES   CELCIUS 

CENTER  TOP  RIGHT  LEFT  BOTTOM  BACK 

CHIPNOI:    28.6rE»00      27.35E*no      27. t3E+00      27.t6E+00      25.93E*00  29.3IE*00 

POWER    (WATTS):       n.S7£-01 

CHIP   N02:    29.73E*00      27.0SE*00      27.l3E*nO      2S.72E*00      27.11E*00  3I.'53E*00 

POWER    (WATTS):      n.53E-ai 

CHIP  N03:  28.47E*00   27.t7E*0Q  27.49E*00   25.82E*00   27. S9E*00    30.6SE*00 
POWER  (WATTS):   n.S4E-01 

CHIPNOI:  28.47E*00  27. 04E*00  27. t^E^OO  26.18E*00  2S.35E»00    30.i:E*00 
POhER  (WATTS):   n,43E-01 

CHIP  N05:  28.21E*00  26.:9E*00  26.42E*00  2G.30E*00  27.03E*00    30.00E*00 
POWER  (WATTS):   n.^8£-01 

CHIP  NOS:  27.4SE*00  27.90E*00  26.98E+00  27.00E*00  27.7SE+00    33.:8E*00 
POWER  (WATTS):   n.A6E-01 

CHIP  N07:  29.12E-^00  27.77E»00  27.03E>00  2S.97E*00  27.25£>00    30.S5E*00 

POWER  (WATTS):   n.SOE-01 

CHIP  N08:  28.5bE*00  27.4SE*00  27.27E*00  27.40E+00  27.18E*00    33.20E*00 
POWER  (WATTS) :   11  .65E-01 

CHIP  N09:  28.I3E*00  27.32E*00  26.S0E*00  27.12E*00  2G.72E*00    32.07E*00 
POWER  (WATTS):   n.55E-01 

HEAT  EXCHANGERS  TEMPERATURES:     RIGHT      LEFT 

BOTTOM:   •  99.44E-01   99.77E-01 

TOP:  I0.05£*00  96.S0E-0t 

BACK  PLANE  TEMPERATURES  ARE  : 

T(55):  19.74E*00 

T(56):  21.13E*00 

1(57):  21.79E*00 

T(72):  21.A2E*00 

T(73):  21.1^E*00 

T(74):  21  .'^(lE^OO  - 

T(75):  21.£3E*00 

T(75):  20.93E^00 

T(77):  20.32E*00 
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TABLE    4 
TEMPERATURE    DATA    FOR    INPUT    POWER   Q^^^    =    1.5    W 

THESE   RESULTS   riRE    STORED    IN   FILE:      D80992015 

AMBIEflT    TEHP   UflS:         23   C 
VULF^'ETER    PE^^DING   UAS :       4.<^UA2    •bl 
BATH    lEh*?    'hi AS'  'U    C 

THIS   RUN   WAb   RECORDED   WITH   FLOW   VIZ 

ALL    TEMPERATURES   ARE    IN   DEGREES   CELCIUS 

CENTER  TOP  RIGHT  LEFT  BOTTOM  BACK 

CHIPNOI:    32.8IE*00      31.19E*00      31.03E*00      3I.08E*00      30.78E*00  34.'^7E*00 

POWER    (WATTS):       U.77E-01 

CHIP   N02:    32.92E*00      30.93E*00      31.02E*00      30.42E*00      31.07E*00  36.44E*00 

POWER    (WATTS):      14.77E-01 

CHIP  N03:  32.36E*00  30.95E*00  31.32E*00  30.50E*00  31.GSE+00    35.33E*00 
POWER  (WATTS):   U.72E-01 

CHIP  N04:  32.S0E*O0  30.86E*00  31.04E>00  29.86E*00  30.79E*00    34.8GE*00 

POWER  (WATTS):   M.58E-01 

CHIP  M05:  32.28E*00  29.79E*00  30.12E*00  29.99E*00  30.81E*00    34,65E*00 
POWER  (WATTS):   14.64E-01 

CHIP  NOG:  30.86E*00  31.77E*00  30.66E*00  30.72E*00  31.43E*00    39.92E*00 

POWER  (WATTS):   14.o2E-01 

CHIP  N07:  33.43E*00  3!.83E*-00  30.89E+00  30.81E*00  31.12E  +  00    35.45E*00 
POWER  (WATTS):   14.80E-01 

CHIP  N08:  32.82E+00  31.45E+00  3l.27E*00  31.40E*00  31.22E>00    38.77E>00 

POWER  (WATTS):   14.87E-01 

CHIP  N09:  32.2IE*00  31.07E*00  30.36E>00  30.97E-00  30.30E*00    37.25E*00 
POWER  (WATTS):   14.74E-01 

HEAT  EXCHANGERS  TEMPERATURES:     RIGHT      LEFT 

BOTTOM:  99.49E-01   99.G9E-01 

TOP:     ■  99.b9E-0l   95.79E-0I 

BACK  PLANE  TEMPERATURES  ARE  : 


T(55): 

-)  t 

.87E*00 

T(5G): 

23' 

,45E^00 

T(57): 

24, 

,33E-00 

T(72): 

23, 

.84E>00 

T(73): 

23. 

.52E*00 

T(74): 

23, 

.83£*00 

T(7S): 

24. 

,  10E*00 

T(76): 

23, 

.25E*00 

T(77): 

23. 

,  10E*00 
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TABLE    5 
TEMPERATURE    DATA    FOR    INPUT    POWER   Q^^    =    1-9    W 


THESE   RESULTS   ARE   STORED    IN   FILE:      D8100M30 


AMBIENT    TEMP  UAS:        23   C 
VOLT'-ETER   READING  WAS:      5.23^-^    *bl 
EATH    TEMP    UAS:  ID    C 

THIS  RUN   WAS   RECORDED   WITH   FLOW   VIZ 

ALL    TEMPERATURES   ARE    IN    DEGREES   CELCIUS 

CENTER  TOP  RIGHT  LEFT  BOTTOM 

CHIPNOI:    37.4SE*00      35.42E*00      3S.24E*00      35.29E*00      34.90E+00 
POWER    (WATTS):      19.10E-01 

CHIP  N02:    37.46E*00     34.63E+00     35.00E*00      34.29E*iJ0     34.99E*00 
POWER    (WATTS):      19.nE-01 

CHIP  N03:    3S.50E*00      34.89E*00      35.33E*00      34.33E*00      35.5IE*00 
POWER    (WATTS):      19.05E-01 

CHIP  N04:    37.22E*00      34.88E*00     35.28E*00      33.81E*00     35.01E*00 
POWER    (WATTS):      18.86E-01 

CHIP   N05:    36.70E*00      33.S5E*00      34.01E*00      33.86E*00      34.97E*00 
POWER    (WATTS):      18.94E-01 

CHIP  NOG:    35.05E-00     36.55E*00     34.59E>00      34.95E*00     3S.73E*00 
POWER   (WATTS):      18.90E-01 

CHIP  N07:    33.33E*00     36.t7E*00     35.t3E*00      34.99E*00     35.34E*00 
POWER    (WATTS):      19.1SE-01 

CHIP  N08:  37.38E*00  35.81E*00  35.42E*00  35.57E*00  35.43E>00 
POWER  (WATTS):   19.23E-01 

CHIP  N09:  35.47E+00   35.22E*00  34.36E-00   35.10E*00   34.40E*00 
POWER  (WATTS):   19.07E-01 


HEAT  EXCHANGERS  TEMPERATURES:     RIGHT      LEFT 

BOTTOM:  99.34E-01   10.02E*00 

TOP:      .  95.29E-0I   92.^7E-01 


BACK  PLANE  TEMPERATURES  ARE  : 


T(55): 

23. 

,54E*00 

T(56): 

25, 

.35E>00 

T(57): 

26. 

.SSE-'-OO 

T(72): 

25, 

.81E*00 

T(73): 

25. 

,3eE>00 

T(74): 

25, 

.72E*00 

T(75): 

25. 

,29E*00 

T(76): 

25, 

.04E+00 

T(77):  24.85E*00 


BACK 

39 

.63E- 

^00 

42 

.43E' 

►00 

40 

.46E' 

►00 

40 

.20£< 

►00 

39 

.a3E^ 

^00 

48, 

.38E^ 

^00 

40, 

,7gE^ 

>00 

45. 

.OSEi 

■00 

43. 

,13E* 

■00 
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TABLE  6 
TEMPERATURE  DATA  FOR  INPUT  POWER  Qj^j^  =  2.3  W 

THESE  RESULTS  ARE  STORED  IN  FILE:   D8100I740 


CHIP  NOt: 
POWER 

CHIP  N02: 
POWER 

CHIP  N03: 
POWER 

CHIP  N04: 
POWER 

CHIP  N05: 
POWER 

CHIP  NOS: 
POWER 

CHIP  N07: 
POWER 

CHIP  N08: 
POWER 

CHIP  N09: 
POWER 


AMBIENT  TE*'?  WAS:    23  C 
VOLTMETER  READING  WAS:   5.8093  *6 I 
BATH  TEMP  WAS:    10  C 
THIS  RUN  WAS  RECORDED  WITH  FLOW  VIZ 

ALL  TEMPERATURES  ARE  IN  DEGREES  CELCIUS 

CENTER     TOP        RIGHT      LEFT 

42.48E*00  39.94E*00  39.75E*00  39.90E*00 
(WATTS):   23.07E-01 

42.44E*00  39.2GE*00  39.57E*00  38.83E>00 
(WATTS):   23.08E-01 

41.90E*00  39.50E>00  40.09E>00  38.99E*00 
(WATTS):   23.00E-01 

42.09E-00  39.2^E*00  29.82E*00  38.09E*00 
(WATTS):   22.78E-01 

4t.49E*00  37.80E*00  38.35E>00  38.t7E*00 
(WATTS):   22.88E-01 

39.G7E*00  41.33E>00  38.98E*00  39.52E>00 
(WATTS):  22.82E-01 

43.53E*00  40.80E*00  39.7tE>00  39.58E*00 
(WATTS):  23.12E-01 

42.48E  +  00  40.3GE  +  00  40.09E-^00  40.20E*00 
(WATTS):  23.22E-01 

41.70E*00  39.86E*00  38.85E+00  39.61E>00 
(WATTS):   23.03E-01 


BOTTOM 
39.39E*00 

39.S1E*00 

40.59E*00 

39.48E*00 

39.60E*00 

40.48E*00 

39.89E-»00 

39.S7E*00 

38.95E*00 


BACK 
45.05E+00 

48.53E*00 

46.17E*00 

45.71E*00 

45.28E*00 

56.04E>00 

46.53E*00 

5t.52E*00 

49.58E*00 


HEAT   EXCHANGERS    TEMPERATURES; 
BOTTOM:  . 
TOP: 


RIGHT  LEFT 

10.00E*00      10.04E+00 
10.ME+00      10.06E*00 


BACK  PLANE  TEMPERATURES  ARE 


T(55): 

25 

S9E+00 

T(B6) 

:   27 

.76E>00 

T(S7): 

29 

42E*00 

1(72) 

:  28 

.37E»00 

T(73); 

27 

87E>00 

1(74) 

:  28 

.23E+09 

T(75): 

29 

04E*00 

T(76) 

•  27 

41E*00 

T(77): 

27. 

I9E*00 
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TABLE    7 
TEMPERATURE    DATA    FOR    INPUT    POWER    Qi^    =    2 . 7    W 

THESE   RESULTS   ARE   STORED    IN   FILE:      081011130 

AMBIENT    TE.^P  WAS:        23  C 
VOLTMETER   READING   UAS:      5.3^4    ♦SI 
BATH    '•EMP    WAS:  10    C 

THIS   RUN   WAS   RECORDED   WITH   FLOW   VIZ 

ALL    TEMPERATURES   ARE    IN   DEGREES    CELCIUS 

CENTER  TOP  RIGHT  LEFT  BOTTOM 

CHIP   NOl:    46.68E+00      43.58E*00      43.37E+00      43.S9E*00      43.06E*00 
POWER    (WATTS):      27.6GE-01 

CHIP   N02:    46.2EE*00      43.02E*00      43.ME*00      42.04E-00      42.86E*00 
POWER    (WATTS):      27. 51E-01 

CHIP   N03:    45.70E+00      43.18E*00      43.85E*00      42.33E*00      44.32E*00 
POWER    (WATTS):      27.S2E-01 

CHIP   N04:    a6.22E-00     42.S8E+00      43.55E*00      41.5SE*00     43.20E*00 
POWER    (WATTS):      27.26E-01 

CHIP   N05:    45.45E*00     40.3lE*00      41.S7E*00      41.55E*00      43.18E+00 
POWER    (WATTS):      27.37E-01 

CHIP  NOS:    43.55E»00     45.08E-^00      4t.98E*00      43.23E>00      44.00E>00 
POWER    (WATTS):      27.29E-01 

CHIP  N07:    47.S6E*00      44.69E*00     43.43E*00      43.20E*00     43.5GE>00 

POWER    (WATTS):      27.55E-01 

CHIP  NOS:  4S.05E*00   43.92E*00   43.56E+00   43.77E*00   43.34E*00 
POWER  (WATTS):   27.78E-01 

CHIP  N09:  44.78E*00   43.22E*00   41.99E>00   42.90E-^00   41.90E  +  Q0 
POWER  (WATTS):   27.55£-01 

HEAT  EXCHANGERS  TEMPERATURES:      RIGHT      LEFT 

BOTTOM:  •  99-.27E-01   99.49E-01 

TOP:  lO.OOE^OO  98.46E-0t 

BACK   PLANE    TEMPERATURES   ARE    : 

T(55):  2S.87E>00 

T(55):  29.09E*00 

T(57):  3J.13E*00 

T(72):  29.93E>00 

1(73):  23.33E*00 

T(74):  29.77E*0&' 

T(75):  30.78E*00 

T(75):  28.32E*00 

1(77):  28.57E*00 


j 

BACK 

49 

.3IE- 

►00 

S3 

.82E- 

►00 

50 

.88E^ 

►  00 

50 

.62E^ 

►00 

49 

.94E^ 

►00 

63, 

.OIEh 

'00 

51, 

.81E^ 

'00 

57, 

,24E^ 

•00 

55. 

,35E^ 

■00 
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TABLE    8 
TEMPERATURE    DATA    FOR    INPUT    POWER   Q^^   =    3.1    W 

THESE   RESULTS   ARE   STORED    IN  FILE:      D8l0nSI5 


AMBIENT    TE^*P   UAS:        23   C 
VCLi.^ETER    READING   UAS:      6.8013    «b1 
BATH    TEf^P   WAS:         10    C 
THIS   RUN   UAS   RECORDED   WITH   FLOU   VIZ 

ALL  TEMPERATURES  ARE  IN  DEGREES  CELCIUS 

CENTER     TOP         RIGHT      LEFT 

CHIP  NOt  :  50.27E*00   46.37E*00   46.S7E*00   4S.36E*00 
POUER  (WATTS):   31.58E-01 

CHIP  N02:  49.47E*00   45.72E*00   4S.94E>00   44.88E»00 

POWER  (UATTS):   31.60E-01 

CHIP  NG3:  48.49E*00   45.7SE*00   46.70E*00   44.97E*00 
POWER  (WATTS):   31.50E-01 

CHIP  N04:  49.72E*00   46.0GE+00   46.89E*00   44.65E*00 

POWER  (WATTS):   31.19E-01 

CHIP  N05:  48.67E*00   43.73E*00   44.88E*00   44.50E>00 
POWER  (WATTS):   31.32E-01 

CHIP  N06:  45.42E*00   48.2SE>00   44.90E*00   46.17E*00 
POWER  (WATTS):   31.23E-01 

CHIP  N07:  5t.67E*00   48.20E>00  46.82E*4)0   46.5IE*00 
POWER  (WATTS):   31.66E-01 

CHIP  N08:  49.44E*00   46.99E*00  46.61E^-00   46.80E*00 
POWER  (WATTS) :   31.80E-01 

CHIP  N09:  48.20E*00   46.J9E+00   44.97E>00   46.00E*00 
POWER  (WATTS):   31.S3E-01 


BOTTOM 

BACK 

46.34E*00 

54.20E*00 

45.65E*00 

58.34E*00 

47.04E*00 

S4.80E*00 

4G.51E*00 

55.03E*0a 

46.^5E*00 

54.10E*00 

47.08E*00 

G8.15E*00 

46.95E*00 

56.86E+00 

46.32E*00 

62.31E*00 

45.25E*00 

60.81E*00 

HEAT  EXCHANGERS  TEMPERATURES: 
BOTTOM: 

TOP: 


RIGHT  LEFT 

10.02E«-00      10.05E*00 
10.09E*00     99.64E-01 


BACK  PLANE   TEMPERATURES   ARE 


T(55) 
T(56) 
T(57) 
T(72) 
T(73) 
T(74) 
T(75) 
T(76) 
T(77) 


28.58E*00 
30.93E*-00 
33.10E*00 
31 .54E+00 
3!  .01E*00 
31 .47E*0a 
32.56E*00 
30.4tiE*0a 
30.21E*00 
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TABLE  9 
REDUCED  DATA  FOR  INPUT  POWER  Qi^  =  0.3  W 


THE  RAW  Emt  DATA  ARE  FROM  THE  FILE:    D80700925 
THE  POWER  SETTING  PER  CHIP  WAS:    0.3  WA^TS 

CHIP      QNET(W)         Tavq-Ts  Nu  "/.UNC  IN  Nu 

1  27.38E-02       56.03E-01       10.61E+00       20.55E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:   104.^3E-02 

•/.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :179.04E-02 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   n0.76E-02 

•/  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   100.88E-02 

2  27.38E-02       57.G0E-01       10.33E+00       20.14E-0I 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-S  IS:   107.4^E-02 

'/.   UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :174.32E-02 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   n0.98E-02 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   t00.88E-02 

3  27.29E-02       57.70E-01       10.28E+00       20.12E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   107.63E-02 

'L   UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :174.03E-02 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   nO.G3E-02 

'L   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   100.88E-02 

4  27.03E-02      55.15E-01       I0.65E+00      20.32E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   102.58E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :182.08E-02 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   109.23E-02 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   100.89E-02 

5  27.14E-02      56.58E-01       10.42E+00      20.42E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *   E-6  IS:   105.'^0E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :1 77. 49E-02 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   109.8GE-02 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   100.89E-02 

6  27.nE-02       B8.1GE-01       10.13E+00      20.00E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *■  E-6  IS:   108.54E-02 

'/.   UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :172.GGE-02 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   109.94E-02 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   100.89E-02 

7  27.45E-02      56.19E-01       10.62E+00      20.52E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:   104.64E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :178.71E-02 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   11!. 08E-02 

'L   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   100.88E-02 

8  27.55E-02  •     58.00E-0l'       10.33E+00       20.04E-0I 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   1Q8.23E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :173.13E-02 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   111.75E-02 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   100.87E-02 

9  27.33E-02      ■57.25E-01       10.37E+00      20.23E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   106.74E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :175.40E-02 

FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   110.70E-02 

'L   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   100.88E-02 


51 


TABLE  10 
REDUCED  DATA  FOR  INPUT  POWER  Qi^i  =  0.7  W 


THE  RAU  Enf  DATA  ARE  FROM  THE  FILE:    D80701535 
THE  POWER  SETTING  PER  CHIP  UAS:    0.7  WATTS 

CHIP      QNET(W)         Tavg-Ts  Nu  7.UNC  IN  Nu 

1  70.30E-02       n.S1E  +  00       13.32E  +  00       13.35E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-5  IS:   228.83E-02 

'L   UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :872.27E-03 

FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   304.72E-02 

'L   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.05E-02 

2  70.32E-02       n.51E+00       13.32E  +  00       13.35E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-S  IS:   228.84E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :872.24E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   30^.80E-02 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   10t.05E-02 

3  70.09E-02       n.51E  +  00       13.28E  +  00       13.35E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-S  IS:   228.78E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :872.A^E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   303.80E-02 

'L   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.06E-02 

4  69.4tE-02       11.36E+00       13.32E+00       13.43E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-G  IS:   225.48E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :883.80E-03 

FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   300.34E-02 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   10t.07E-02 

5  69.70E-02       11.23E+00       13.53E+00       13.49E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   222.G4E-02 

'L   UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :893.83E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   30t.13E-02 

1   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.06E-02 

6  G9.61E-02       11.34E+00       13.38E+00       t3.44E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   22^.97E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :885.5GE-03 

FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   301.12E-02 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.07E-02 

7  70.49E-02       11.57E+00       13.28E+00       13.32E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   230.19E-02 

'L   UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :867.G9E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   305.75E-02 

•/  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.05E-02 

8  70.78E-02  •     n.62E-t-00''      13.28E+00       13.30E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   231.24E-02 

7  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :8G4.18E-03 

FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   307.15E-02 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.05E-02 

9  70.17E-02       41.34E*00       13.49E+00       13.43E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   225.08E-02 

'L   UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :885.19E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   303.58E-02 

1   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.06E-02 
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TABLE  11 
REDUCED  DATA  FOR  INPUT  POWER  Q^^    =  1.1  W 


THE  RPW  Ent-  DATA  ARE  FROh  THE  FILE:    080711035 
THE  FOUER  SETTING  PER  CHIP  WAS:    1.1  WATTS 

CHIP      ONET(W)         Tavg-Ts  Nu  '/.UNO  IN  Nu 

1  n.45E-01       16.67E^00       15.01E  +  00       11.77E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   343.22E-02 

2   UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :602.83E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   515.31E-02 

•/  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.06E-02 

2  11.4GE-0I       15.59E+00       15.08E+00       n.78E-0I 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   3^1.44E-02 

•/.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :G05.54E-03 

FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   51S.06E-02 

•/.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.08E-02 

3  11.42E-01       16.56E  +  00       IS-OGE  +  OO       n.79E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   340.72E-02 

•/.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :605.64E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   513.2tE-02 

X  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.06E-02 

4  n.31E-01       15.41E+00       15.05E  +  00       n.82E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   337.05E-02 

•/.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :6J2.23E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   507.35E-02 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.07E-02 

5  11.35E-01       15.05E+00       15.45E+00       n.89E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-G  IS:   328.29E-02 

2   UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :G26.12E-03 

FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   507.35E-02 

2   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.07E-02 

S       n.34E-01       15.92E+00       15.55E+00       n.92E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   32S.25E-02 
7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :631.nE-03 
FLUX  BASED  RAYLEIGH  NUMBER  •  E-7  IS:   505.88E-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.07E-02 

7  11.48E-01       16.81E+00       14.93E+00       11.74E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   346.G8E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :597,78E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   517.48E-02 

•/  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.06E-02 

8  n.53E-01  ■     16.76E+00       15.03E  +  00       n.75E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   345.42E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :599.62E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   B19.33E-02 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   t01.05E-02 

9  11.43E-01     ^  1G.28E^-00       15.34E  +  00       n.84E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:   333.97E-02 

•/  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :617.04E-03 

FLUX  BASED  RAYLEIGH  NUMBER  •»  E-7  IS:   512.18E-02 

1   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.06E-02 
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TABLE  12 
REDUCED  DATA  FOR  INPUT  POWER  Q^j^  =  1.5  W 


THE  RAU  Emf  DAT^  AFE  FROM  THE  FILE:    D80720730 
THE  POWER  SETTING  PER  CHIP  WAS:    1.5  WATTS 

CHIP      QNET(W)         Tavg-Ts  Nu  7.UNC  IN  Nu 

1  M.59E-01       19.05E^00       16.76E  +  00       n.41E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   401.43E-02 

1   UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :527.G7E-03 

FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   S72.7!E-02 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   10*.12E-02 

2  I4.59E-01       I9.08E+00       t6.7^E+00       11.40E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-b  IS:   402.25E-02 

•J.   UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :52G.77E-03 

FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   S73.17E-02 

•/.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.t2E-02 

3  14.SAE-01       19.35E*00       16.A5E+00       n.37E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   409.23E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :519.31E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   G73.03E-02 

•/  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.12E-02 

4  U.40E-01       t8.73E  +  00       1S.83E  +  00       n.45E-0I 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   393.35E-02 

•/.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :536.57E-03 

FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   661.82E-02 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.13E-02 

5  14.46E-01       18.39E+00       17.20E+00       11.49E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-5  IS:   385.00E-02 

X   UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :5A6.38E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   662.05E-02 

'L   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.13E-02 

S        I4.44E-01       18.59E+00       1G.99E+00       11.47E-0t 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-5  IS:   389.8GE-02 
7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :540.G8E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   662.4GE-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   10T.13E-02 

7  14.62E-01       19.27E+00       15.G0E+00       11.38E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   407.23E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :521.A2E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   67G.12E-02 

•L   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   IOI.nE-02 

8  14.68E-01  ■     19.48E>00       1G.50E+00       n.35E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-5  IS:   412.33E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :5t6.06E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   G80.^7E-02 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.11E-02 

9  14.56E-01      -ig.OgE+OO       1G.G9E+00       11.A0E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-5  IS:   402.G0E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :526.39E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   671.79E-02 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.12E-02 
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TABLE  13 
REDUCED  DATA  FOR  INPUT  POWER  Q^^   =  1.9  W 


THE  RAW  E.Ttt  DATA  ARE  FROM  THE  ^ILE:    D80721810 
THE  POUER  SETTING  PER  CHIP  UAS:    1.9  WATTS 

CHIP      OMET(W)         Tav/g-Ts  Nu  XUNC  IN  Nu 

1  18.99E-01       24.92E+00       16.75E+00       10.88E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   565.77E-02 

y.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :403.27E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ■*   E-7  IS:   947.59E-02 

2   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.07E-02 

2  19.00E-01       2'!l.72E  +  00       15.89E  +  00       10.89E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   5G0.03E-02 

v.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :406.56E-03 

FLUX  BASED  RAYLEIGH  NUMBER  •  E-7  IS:   945.98E-C2 

•/.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.07E-02 

3  18.9AE-01       2/J.58E  +  00       16.86E  +  00       10.90E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   559.04E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :407.13E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   942.52E-02 

1   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.08E-02 

A  18.75E-01       24.53E+00       15.80E+00       10.91E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   55^.73E-02 
•/  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :'^09.65E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   931.79E-02 
y.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.09E-02 

5       18.83E-01       23.88E+00       17.32E+00       10.95E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-5  IS:   536.40E-02 
7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :420.80E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ■»  E-7  IS:   928.99E-02 
X   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.08E-02 

e       18.80E-01       23.83E+00       17.32E*00       10.95E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   535.00E-02 
7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :421.68E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   926.8AE-02 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   10t.09E-02 

7  19.04E-01       2B.23E+00       16.59E+00       10.8SE-01 
TEMPERATURE-  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   574.77E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :398.25E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   953.30E-02 

•/  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.07E-02 

8  19.t2E-01  ■     25.20E+00       16.58E+00       10.8GE-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-S  IS:   573.88E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :398.74E-03 

FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   957.00E-02 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   10t.07E-02 

9  18.96E-01      •2^.47E+00       17.02E+00       10.91E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-5  IS:   552.91E-02 

Z   UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :410.73E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   941.24E-02 

X   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.08E-02 
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TABLE  14 
REDUCED  DATA  FOR  INPUT  POWER  Qj^n  =  2.3  W 


THE  RAW  Enf  DATA  ARE  FROM  IHE  FILE:    D80731030 
THE  POWER  SETTING  PER  CHIP  WAS:    2.3  WATTS 

CHIP      QNET(W)         Tavg-Ts  Nu  7.UNC  IN  Nu 

1  22.85E-01       28.31E+00       17.76E+00       10.71E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   658.S3E-02 

•/  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :3S5.02E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   :i6.98E-01 

'L   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.07E-02 

2  22.86E-01       27.97E+00       17.98E+00       t0.73E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   6A8.A0E-02 

•/.   UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :359.32E-03 

FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   IIS.SOE-OI 

•/.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   tO!.07E-02 

3  22.78E-01       27.97E+00       17.93E+00       10.73E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  -*  E-S  IS:   648.38E-02 

'L   UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :3S9.34E-03 

FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   nS.22E-01 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.07E-02 

4  22.56E-01       27.37E+00       17.SIE+00       10.73E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-5  IS:   G45.4GE-02 

•/.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :3S0.59E-03 

FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   1M.37E-01 

v.   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.08E-02 

5  22.G5E-01       27.10E+0n       18.38E+00       10.77E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-S  IS:   G22.74E-02 

•/  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :370.80E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   1M.48E-01 

'L   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.08E-02 

G       22.60E-01       27.!5E*00       18.31E+00       10.76E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-G  IS:   G24.2GE-02 
7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :370.10E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   114.28E-01 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.08E-02 

7  22.90E-01       28.72E+00       17.55E+00       10.70E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-G  IS:   G71.04E-02 

•/  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :349.89E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   n7.80E-01 

•/  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.07E-02 

8  23.00E-01  ■     28.59E+00       17.71E+00       10.70E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-G  IS:   6G7.05E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :351.51E-03 

FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   118.13E-01 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.0GE-02 

9  22.80E-01      ■27.G8E+00       18.13E+00       10.74E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-G  IS:  639.69E-02 

•/  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :3G3.12E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-/  IS:   I15.97E-0I 

'L   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.07E-02 
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TABLE  15 
REDUCED  DATA  FOR  INPUT  POWER  Q^^    =  2.7  W 


THE  RAW  Emt  DATA  ARF  FROM  THE  ^ILE:    DSn731755 
THE  POWER  SETTING  PER  CHIP  WAS:    2.7  uATTS 

CHIP      QNET(W)         Tavg-Ts  Nu  7.UNC  IN  Nu 

1  27.37E-01       33.28E+00       18.16E+00       1C.55E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  >  E-6  IS:   313.30E-02 

•/.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :302.08E-03 

FLUX  BASED  RAYLEIGH  NUMBER  *  £-7  IS:   147.70E-01 

X   UNCERTAINTY  IN  FLUX. BASED  RAYLEIGH  NUMBER  IS:   101.05E-02 

2  27.39E-01       32.77E+00       18.44E+00       10.56E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-S  IS:   79G.95E-02 

•/  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :306.72E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   146.98E-01 

•/.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.05E-02 

3  27.30E-01       32.53E+00       18.52E+00       10.57E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:   7S9.23E-02 

•/.   UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :308.97E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   UG.t3E-0t 

2   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.05E-02 

A  27.03E-01       32.75E+00       18.22E*00       10.56E-01 

TEMPERATURE  BASED  RAYLEIGH  NUMBER  *  E-6  IS:   796.28E-02 
7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :306.91E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   145.04E-01 
7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.06E-02 

5  27.15E-01       31.73E+00       18.87E+00       10.59E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *   E-5  IS:   7S3.51E-02 

1  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :315.79E-03 
FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   144.06E-01 

•/.   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101,06E-02 

6  27.07E-01       31.79E+00       18.78E+00       10.59E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-5  IS:   765.48E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :316.18E-03 

FLUX  BASED  RAYLEIGH  NUMBER  *  E-7  IS:   1A3.75E-01 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   I01.06E-02 

7  27.44E-01       33.82E-^00       17.92E+00       10.53E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  *   E-S  IS:   83t.ME-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :297.22E-03 
FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   148.91E-01 

2  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.05E-02 

8  27.56E-01  '     33.41E-^00       18.21E+00       10.54E-0t 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-6  IS:   817.5GE-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :300.90E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   1A8.90E-01 

7.  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.04E-02 

9  27.32E-01      -32.11E+00       18.77E+00       10.58E-01 
TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-5  IS:   775.54E-02 

7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  IS  :313.07E-03 

FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   MB.59E-01 

1   UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   101.05E-02 
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1.   Graphical  Representation  of  Reduced  Data 

The  experimental  surface  temperature  data  were 
plotted  as  a  dimensionless  Nusselt  Number  vs.  both  a 
temperature  based  and  a  flux  based  Rayleigh  Number,  Figures 
4.15  through  4.32. 

Representative  error  bounds  were  plotted  in  Figures 
4.15  and  4.16  for  input  powers  of  1.5  and  1.9  watts.  All 
points,  with  the  exception  of  1.5  watts  were  well  within  the 
error  bounds.  The  point  representing  1.5  watts  was  either 
too  high  on  the  Nusselt  Number  axis  or  too  low  on  the 
Rayleigh  number  axis. 

The  factors  affecting  the  Nusselt  Number  were  h,  L  & 
K.  Of  these  three,  only  h  was  calculated  with  uncertainty 
in  this  experiment.  The  factors  affecting  the  Rayleigh  Num- 
bers were  a,  L,  AT,  Qnet'  3'  v  ^rid  K.  Of  these  only  AT  and 
Qnet  ^^^  uncertainties. 

Examination  of  Tables  4  and  12  resulted  in  no  data 
which  appeared  to  be  extreme  in  either  the  raw  or  reduced 
data  sets.  One  possible  explanation  for  this  point  is  the 
changing  flow  patterns  around  this  power  input.  This  may  be 
a  transition  region  from  somewhat  orderly  flow  to  completely 
chaotic  flow.  It  is  recommended  that  further  investigation 
of  this  region  of  power  settings  be  conducted. 

The  equations  for  the  best  fit  lines  are  very  close 
to  each  other  suggesting  that  an  overall  correlation  that 
would  produce  good  approximations  is  possible. 
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Figure  4.15   Temperature  Based  Rayleigh  Numbers,  Chip  1 
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Figure  4.16   Flux  Based  Rayleigh  Numbers,  Chip  1 
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CHIP  NUMBER  2 


TEMP  BASED  RAYLEICH  NUMBER 


Figure  4.17   Temperature  Based  Rayleigh  Numbers,  Chip  2 


CHIP  NUMBER  2 


20 

j      i 

i    i 

i  1 

; 

I 

■■]-Ti  \ 

y^  I 

19 

1 

\    \    1 

y^ 

i 

y^^ 

IB 

1 

: 

I 

; 

y^ 

17 

- 

1 

i 
oi 

y^ 

1 

UI6 

m 

^ 

: 

^ 

i 

* 

5,. 

10 

I/) 

D 
ZI3 

. 

y 

^ 

i 

i 

12 

y 

r^ 

\l 

j 

II 

^o 

1  " 

\ 

-to 

o.ttto 

Nw  0.2606  Ro 

"l 

0' 

' 

1 

\i 

FLUX  BASED  RAYLEICH  NUMBER 


Figure  4.18   Flux  Based  Rayleigh  Numbers,  Chip  2 
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Figure  4.19   Temperature  Based  Rayleigh  Numbers,  Chip  3 
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Figure  4.20   Flux  Based  Rayleigh  Numbers,  Chip  3 
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Figure  4.21   Temperature  Based  Rayleigh  Numbers,  Chip  4 
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Figure  4.22   Flux  Based  Rayleigh  Numbers,  Chip  4 
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Figure  4.2  3   Temperature  Based  Rayleigh  Numbers,  Chip  5 
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Figure  4.24   Flux  Based  Rayleigh  Numbers,  Chip  5 
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Figure  4.25   Temperature  Based  Rayleigh  Numbers,  Chip  6 
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Figure  4.26   Flux  Based  Rayleigh  Numbers,  Chip  6 
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Figure  4.27   Temperature  Based  Rayleigh  Numbers,  Chip  7 
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Figure  4.2  8   Flux  Based  Rayleigh  Numbers,  Chip  7 
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Figure  4.29   Temperature  Based  Rayleigh  Numbers,  Chip  8 


CHIP  NUMDER  8 


20 
19 

18 
IT 

UJ  16 

m 
I" 


to 
ZI3 


12 


1 

i 1 1 1 

i       i 
]       I 
i       1 

- 

1            ;         i 

i    i    1 

1         1       ; 
1         ■       ■ 

i 

\"\ 

i       I 

N 

/^o 

■ 

i        i      i 

i     1    1 

I        1      i 

:                  !              ■         > 

!              i           :a/ 

O 

!     ■ 
t 

i 

! 

-^ 

1 

- 

Jr         \        i   \ 

y\      \    \  \ 

:                i            !         i 

\ 

i 

1 

1 

'o 

i                        :                 :             i 

r                                         : 

:                                            <              : 

i                           !         i 

;                        !                 i             i 

:V 

Nu»  0  2953  Ro 

I'o' 

. 

1 

'(/ 

FLUX  DASEU  RAYLEIGH  NUMDER 


Figure  4.30   Flux  Based  Rayleigh  Numbers,  Chip  8 
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Figure  4.31   Temperature  Based  Rayleigh  Numbers,  Chip  9 
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Figure  4.32   Flux  Based  Rayleigh  Numbers,  Chip  9 
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Figure  4.33   Temperature  Variation  for  Input  Power  0.3  W 
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Figure  4.34   Temperature  Variation  for  Input  Power  0.7  W 


69 


C3> 


^   o   o 

(\J    N    CM 


o 

"  0 

o  • 

E 

2  ^ 


N 


rg 

^ 

C3> 

^       D       o 

(\J          f\)          N 

o 

10 

I\l 

ID 
(M 

Figure  4.35  Temperature  Variation  for  Input  Power  1.1  W 
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Figure  4.3  6   Temperature  Variation  for  Input  Power  1.5  W 
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Figure  4.37  Temperature  Variation  for  Input  Power  1.9  W 
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Figure  4.38   Temperature  Variation  for  Input  Power  2.3  W 
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Figure  4.39  Temperature  Variation  for  Input  Power  2.7  W 
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Figure  4.40   Temperature  Variation  for  Input  Power  3.1  W 
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C.   TEMPERATURE  VARIATIONS 

The  temperature  variations  with  time  were  investigated 
and  reported  in  Figures  4.3  3  through  4.40.  These  results 
were  motivated  by  the  studies  conducted  by  Liu,  Yang  and 
Kelleher  [Liu  et  al.,  1987,  pp.  486-97].  Three  thermocouple 
channels,  numbers  1,  13  and  32  were  scanned  individually  at 
a  rate  of  3  times  per  second  for  a  period  of  200  seconds. 
These  quantitative  results  were  then  plotted.  The  results 
clearly  demonstrate  a  multiple  frequency  content  signal  for 
most  of  the  power  levels  investigated  in  this  experiment. 

An  observation  of  this  experiment  was  that  the  output  of 
thermocouple  #1  remained  steady  with  time  throughout  the 
increase  in  power  levels.  The  temperature  level  for  this 
thermocouple  was  also  generally  higher  than  for  the  other 
two  thermocouples  sampled. 

A  possible  explanation  for  these  results  is  that  the 
cold,  stagnant  fluid  layer  created  by  the  bottom  heat 
exchanger  prevented  the  penetration  of  warmer  fluid  below 
the  bottom  row  of  components.  This  is  also  supported  by 
flow  visualization. 

One  exception  is  at  3.1  watts.  As  seen  in  Figure  4.40, 
there  is  a  large  fluctuation  of  temperature  with  time.  A 
possible  explanation  for  this  fluctuation  may  be  that  the 
increased  energy  acquired  by  the  fluid  at  this  power  setting 
enabled  the  flow  to  force  itself  into  the  stagnant  region  at 
the  bottom  of  the  chamber. 
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V.   RECOMMENDATIONS 

There  is  a  great  deal  more  experimentation ,  to  be 
accomplished  with  this  component  configuration.  Three 
possible  variations  to  the  chamber  are  : 

-  Placing  a  spacer  in  the  chamber  to  reduce  the  chamber 
width.  This  would  act  in  the  same  manner  as  a  shrouding 
wall  and  would  affect  the  heat  transfer  characteristics. 

-  Setting  the  bottom  heat  exchanger  to  a  temperature  other 
than  10°C.  The  energy  expended  to  maintain  the  bottom 
heat  exchanger  at  10°C  is  considerable.  If  reasonable 
components  temperatures  can  be  realized  by  maintaining 
the  bottom  heat  exchanger  at  a  higher  temperature  level 
then  the  energy  consumption  of  the  system  would  be 
reduced. 

-  Powering  different  configurations  other  than  the  entire 
array.  This  could  possibly  offer  some  explanation  of  the 
flow  patterns  observed  in  the  chamber. 

Other  possible  improvements  to  software  and  hardware 


are: 


Producing  a  plotting  program  that  would  read  the  data 
directly  from  the  data  acquisition  unit  for  plots  of 
Nusselt  vs.  Rayleigh  number. 

Implementation  of  measurement  devices  that  have  greater 
sensitivity  to  minor  fluctuations  in  temperature 
variation. 

Development  of  a  program  that  will  perform  frequency 
analysis  on  the  data  recorded  by  the  program  Fastscan. 
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APPENDIX  A 
SAMPLE  CALCULATIONS 

Throughout  this  appendix,  the  sample  calculations  will 
be  based  on  component  number  one  with  a  power  input  of  0.71 
watts. 
1.   Determination  of  Input  Power 


^•r'h }-  h  \-  \-l-i- 


rnTTT 

V.   V^   V    V   V   V   V   V   V 

I    i  1 1  1  i  11  ii 


Figure  A. 1   Electrical  Network  of  Power  Input 


From  Equation  (3.4): 


Qin  =  ((3.218  -  0.5278) (0.5278) )/2 


78 


2.  Estimation  of  Conduction  Loss 

Equation  3.3  requires  T^,  T^  and  Rcon*  "^s  ^^^  "^b  ^©re 
read  directly  from  thermocouples  placed  immediately  behind 
the  heaters  and  on  the  back  of  the  board,  respectively, 
^con  ^^^  calculated  by  using  the  constriction  resistance 
formulas  developed  by  Mikic  and  repeated  in  [Kraus  and  Bar- 
Cohen,  1983,  pp.  76-9].  The  conduction  losses  for  components 

3 ,  5  and  7  were  then  averaged  to  yield  a  total  conduction 
loss. 

3 .   Calculation  of  Nusselt  and  Rayleigh  Numbers 

The  calculations  of  the  convection  coefficient  from 
Equation  (3.5)  was: 

Qnet  =  hA(Tavg  "  T^) 

Qnet 
""   =  ^(^avg  -  ^c) 

.703 


(5.76  X  10"^)  (11.51) 


=  106.04  W/m2  Or 

From  Equation  (3.10): 

T     +  T 
avg    c 

Tfilm  =  2 

21.51  +  10 
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Tfilm  =  15.75OC 

Using  the  relationships  of  [3-M  Manual,  1985,  pp.  9-60],  the 
physical  properties  of  FC-75  were: 

kf  =  (.65  -  7.8947  X  10""^  *  Tfiiin)(0-1) 

=    0.0637    W/in°C 

p      =    1000(1.825    -    0.00246    *    Tfi^m) 
=    182  3.21   Kg/in^ 

Cp    =    4180(. 241111    +    3.7037    X    10""^    *    Tfii^) 
=    1032.23    J/Kg°C 

V  =    10"^(1.4074    -    2.964    x    10"^    *    Tfiij^ 

+  3.8018  X  10""*  *  Tfiiju 
-  2.7308  X  10"^  *  Tfilm 
+    8.1679    X    10"^    *    Tfiijii) 

V  =   1.024   X   10"^  m^/sec 


.00246 

^  1.825    -    .00246    *    T^., 

rilm 


3     =      1.377    X    10 


-3 
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Then  from  Equation  (3.1): 


Nu  =  ^ 


33  106.04  (  .008) 
.0637 


=  13.32 


K 

a  = 


pc 
P 


.0637 


1823.2  *  1032.2 


=  3.38  X  10"^  m^/sec 


Pr  =  ^ 

a 


1.024  X  10  ^ 
3.38  X  lO"^ 


=  30.3 


Grt  -         2 


9.81(1.377  X  10"^) (.008)^ (11.51) 
(1.024  X  10~^)^ 


=  7.592  X  10^ 
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From  Equation  (3.11),  the  temperature  based  Rayleigh  number 
was: 

Rat  =  Grt  *  Pr 
=  2.3  X  10^ 

From  Equation  (3.12),  the  flux  based  Rayleigh  number  was: 


Raf  =  3.06  X  10*7 
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APPENDIX  B 
UNCERTAINTY  CALCULATIONS 

Uncertainty  calculations  were  evaluated  using  the  root 
mean  square  formulation  described  in  [Beckwith,  Buck  and 
Marangoni,  1982,  pp.  261-292].  The  following  calculations 
are  for  components  number  one  at  0.7  watts  (see  Table  4.10): 

TABLE  16 
UNCERTAINTIES  IN  VARIABLES 


Variables 

Uncertaint 

ies 

Reference 

L 

10~5  m 

Resolution 

T 

0.1  Oc 

Thermocouple  readings 

Rp 

0.02 

1%  Variation 

Kp 

1.5  % 

[Touloukian,  1986,  p.  970] 

V 

0.001  V 

Resolution 

The  uncertainty  in  the  net  power  added  to  the  fluid  was; 


6Q   .  =  (6Power   +  6Q^^^^  ) 


2,1/2 


where 


6  Power  =  0.01  *  Power 


=  .0071  W 
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SQloss  =  0.055  *  Qioss 
=  .000159  W 


Therefore, 


6Qnet  =  -OO^l 


The  uncertainty  for  the  other  variables  used  were; 


^Tsink  =  «T 


=  0.1  K 


^Tfilm  =    (5T|ink  +    ^avq)  ^^^ 


=    (.0l2    +    0.l2)V2 


=    .1414    K 


SToelt  =    («Tavg  +   «^sink)  ^^^ 
=    (0.l2    +    .012)1/2 


=    .1414    K 


84 


5h 
h 


^^net  2     ^T        1/2 


"net 


Delt 


r  ,.0071  2  _^  ,.1414>  2. 


1/2 


.703 


11.51 


=  .0159 


6Nu 
Nu 


6h_ 
h 


=  .0159 


6Gr^    6T 


Delt 


Gr, 


Delt 

.1414 
11.51 


=  .0123 


~RaT 


5Gr^ 
"GrT 


=  .0123 
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5Ra 
Ra 


f  _ 


6Q 


net 


'net 

.0071 
.703 


=  .0101 
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APPENDIX  C 
SOFTWARE 


1.   Acquire 


10 
20 
30 
40 
A\ 
42 
SO 
60 
70 
80 
90 

too 

110 
120 
130 
140 

150 

160 

170 

180 

190 

200 

210 

211 

220 

230 

240 

2S0 

260 

261 

270 

280 

290 

291 

300 

310 

320 

330 

340 

3S0 

360 

370 

380 

390 

400 

10ft" 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 


'•,10fl";01df  ile$ 


10A*:Newf jle$ 


!  FILE  ACQUIRE 
I 

!  EDITED  BY  LT  T  J  BENEDICT.  USN 
I 

!  READ  FILE  'READ  ME' 

COH  /Co/  0(7) 

DIM  Emf (76).Power(9) ,T(75) 
I 

ICORRELATION  FACTORS  TO  CONVERT  EMF  TO  DEGREES  CELCIUS 
DATA  0.10086091 .25727.9.-767345.8.78025596, 
DATA  -924  7486589. 6. 98E  I  1  . -2.56E  I  3.3.  94E  l'4 
I 

READ  D(«) 

Rp-2.000 

PRINTER  IS  701 

BEEP 

! 

INPUT  'ENTER  THE  INPUT  MODE:  0-SYS,  l-FILE'.Im 
! 

IF  Im-1  THEN 

BEEP 

INPUT  'ENTER  THE  NAME  OF  THE  FILE  TO  BE  READ" .Oldf l le$ 
I 

PRINT  USING  'ISX, "THESE  RESULTS  ARE  STORED  IN  FILE  : 
I 

ELSE 

BEEP 

INPUT  "ENTER  THE  NAME  OF  THE  NEW  FILE" .Newf i le$ 

PRINT  USING  '10X, "THESE  RESULTS  ARE  STORED  IN  FILE:   ' 

END  IF 

PRINT 

PRINT 
j 

INPUT  'ENTER  THE  AMBIANT  TEMP" . TemoS 

PRINT  USING  '15X."'AMBI£NT  lEMP  HAS:    •" , 10A" : Te«p$ 

I 

Input  "Flou  viz?  Y/N",Ans$ 

INPUT    "VOLTMETER  READING*. V$ 

PRINT   USING    'ISX, "VOLTMETER  READING   WAS:      ".lOfl'lVS 

I 

INPUT  "ENTER  THE  BATH  TEMP'.BS 

PRINT  USING  "ISX, ""BATH  TEMP  WAS    :"'.10A':B$ 

I 

IF  Ans$-"Y"  THEN  PRINT  USING  "ISX. ""THIS  RUN  WAS  RECORDED  WITH  FLOW  VIZ'", 

I 

IF  Im-1  THEN  ASSIGN  «File  TO  Oidfile« 

PRINT 

! 

IF  Im-0  THEN 

CREATE  BOAT  Newflle$.S 

ASSIGN  «Flle  TO  Newf lle$ 

END  IF 

I 

!  READ  DATA 

I 

IF  Im-0  THEN 

OUTPUT  709; "AR  AFOO  AL79" 
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S40  OUTPUT  72r:T1  R1  T1  ZO  FLO* 

5S0  ! 

560  FOR  I-O  TO  7G 

S70  OUTPUT  709: "OS" 

580  UAIT  I 

590  ENTER  722:Enif(I) 

600  BEEP 

610  NEXT  I 

620  ! 

630  OUTPUT  »Fj le:Emf (•) 

640  ! 

6S0  ELSE 

650  EMTER  *F  i  le'.Emf  (  •) 

f^70  END  IF 

680  ! 

690  OUTPUT  709:"TD" 

700  ! 

710  FOR  1-0  TO  60 

720  Sum-0 

730  FOR  J-1  TO  7 

740  Sum-Suin*D(  J)»Einf  (I)*  J 

7S0  NEXT  J 

760  T(I)-Sum 

770  NEXT  I 

780  ! 

790  FOR  r-71  TO  76 

800  Sum-0 

810  FOR  J-1  TO  7 

820  Sum-Suin*D(  J)«Emf  (  I)  •  J 

830  NEXT  J 

840  T(I)-Sun 

OSO  NEXT  1 

860  ! 

870  I  POWER  CALCULATIONS 

880  ! 

890  J-1 

900  Volt-Einf(6n 

910  • 

920  FOR  1-62  TO  70 

930  Power(J)-Einf  (I)»(Volt-E«f  <I))/Rp 

940  J-J*I 

9S0  NEXT  I 

960  ! 

970  BEEP 

980  BEEP 

990  PRINT  USING  'MX,  ""ALL  TEMPERATURES  ARE  IN  DEGREES  CELCIUS'"" 

1000  ! 

1010  PRINT 

1020  ! 

1030  PRINT  USING  " 12X, ""CENTER     TOP         RIGHT      LEFT      BOTTOM 

BACK"" 

1040  PRINT 

1050  PRINT  USING  "IX, ""CHIP  N01:  "" .5( ZZ.DDE .2X) .2X.ZZ.DDE" : T (0) . T( 1 ) . T(2) , T( 
.T(4) ,1(5) 

1060  PRINT  USING  "SX, ""POWER  (WATTS):   "" , ZZ. DDE" :Power (1 ) 

1061  PRINT 

1070  PRINT  USING  "IX. ""CHIP  N02:  "" .5( ZZ.DDE .2X) .2X, ZZ.DDE" ; T (6) , T ( 7 ) . T (8) , T( 
,T(10).T(in 

1080  PRINT  USING  "5X, ""POWER  (WATTS):   "", ZZ.DDE" :Power (2) 

1081  PRINT 
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1090  PRINT 
(  IS)  .T(  \f))  , 

1100  PRINT 

1101  PRINT 
nin  PRINT 
(21 ) .T(22) . 

1120  PRINT 

1121  PRINT 
1130  PRINT 
(27)  .T(28) . 
1  MO  PRINT 
11A1  PRINT 
nSO  PRINT 
(33) .T(34). 

n&O  PRINT 

1161  PRINT 

1170  PRINT 
(3?)  ,  TC^O)  . 

1180  PRINT 

1181  PRINT 
Itgo  PRINT 
<  ^S)  .  T(/t6  )  . 

1200  PRINT 

1201  PRINT 
1210  PRINT 
(SI  ).T(S2). 
1220  PRINT 
1230  ! 
12^0  PRINT 
nsp  PRINT 

1270  PRINT 

1280  PRINT 

1290  PRINT 

1300  P»INT 

1310  PRINT 

1320  ! 

1330  PRINT 

1340  PRINT 

13S0  PRINT 

1360  PRINT 

1370  PRINT 

1380  PRINT 

1390  PRINT 

1400  PRINT 

1410  PPINT 

1420  PRINT 

1430  PRINT 

1440  BEEP 

14S0  PRINTER 

1460  • 

1470  ASSIGN  ^File 

1480  END 


USING  -IX." 
T(17) 

USING  "SX." 

USING  'IX.- 
1(23) 

USING  "SX," 

USING  "IX." 
T(29) 

USING  'SX.' 

USING  "IX," 
T(3S) 

USING  'SX.- 

USING  'IX.- 
T(41  ) 

USING  'SX.- 

USING  "IX.- 
T(47) 

USING  -SX.- 

USING  MX,' 
T(53) 

USING  'SX." 


USING 
USING 
USING 


CHIP  N03:  ■•.S(ZZ.DDE.2X).2X,ZZ.DDE';T(I2).T(13).T(14) 
POUER  (WATTS):   ••.ZZ.DDE':Power(3) 

CHIP  N04:  ••.5(ZZ.nDE.2X) .2X.ZZ.DDE':T( I8).T( 19) .1(20) 
POWER  (WATTS):   "  , ZZ . DDE" :Power ( 4 ) 

CHIP  NOS:  ".5(ZZ.DDE.2X).2X.ZZ.DDE":T(24),T(2'S).T(26) 
POWER  (WATTS):   ■•.ZZ.DDE":Power(S) 

CHIP  N06:  ••,S(ZZ.DDE.2X).2X.ZZ.DDE';T(30),T(31).T(32) 
POWER  (WATTS):   "  . ZZ . DDE" :Power (6) 

'CHIP  N07:  •".S(ZZ.DDE.2X).2X.ZZ.DDE":T(36).T(37) .T(38) 
POWER  (WATTS):   "" . ZZ . DDE" :Power ( 7 ) 

CHIP  N08:  ••.5(ZZ.DDE.2X).2X.ZZ.DDE":T(42).T(43).T(44) 
■POWER  (WATTS):   ■•.ZZ.'DDE":Power(8) 

CHIP  N09:  ••.S(ZZ.DDE,2X).2X.ZZ.DDE":T(48).T(49).T(50) 
POWER  (WATTS):   "".ZZ.DDE":Power(9) 


'5X." 

'10X. 

10X, 


HEAT  EXCHANGERS  TEMPERATURES:     RIGHT 
'•B0TT0M:"*.21X.2(ZZ.DDE.2X)":T(S8).T(S7) 

•TOP:"".24X.2(ZZ.DDE.2X)":T(60).T(S9) 


LEFT' 


USING  •SX,""BACK  PLANE  TEMPERATURES  ARE 


USING 
USING 
USING 
USING 
USING 
USING 
USING 
USING 
USING 


IS  I 


10X. 
'10X. 

10X. 
'10X, 
'10X. 
'10X. 
'10X. 
'10X, 
'10X. 


TO  • 


"T(5S) 
'•T(S5) 
'•T(S7) 
••T(72) 
'•T(73) 
'•T(74) 
"T(75) 
•"T{76) 
'"T(77) 


.2X.(ZZ.DDE.2X1" 
.2X.(ZZ.DDE,2X)" 
.2X,(ZZ.DDE.2X)" 
.2X.(ZZ.DDE.2X)" 

.2X.(ZZ.DDE.2X)" 
.2X.(ZZ.DDE,2X)' 
.2X.(ZZ.DDE.2X)" 
.2X,(ZZ.D0E.2X)* 
,2X.(ZZ.DDE.2X)" 


T(S4) 
T(S5) 
T(S6) 
T(71) 

T(72) 
T(73) 
T(74) 
T(7S) 
T(76) 
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2.   CalcDiel 


10 

20 

30 

AO 

SO 

80 

70 

80 

90 

100 

110 

120 

130 

MO 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

2S0 

260 

270 

280 

290 

300 

310 

320 

340 

350 

370 

380 

390 

400 

A10 

420 

430 

440 

450 

450 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

610 

520 


FILE  NrtflE  :  CdlcDiei 

EDITED  BY  :  Lt.  T  BENEDICT.  USN 

THIS  PROGRAM  ANALYSES  THE  DATA  READ  FROM 
A  DATA  FILE  DESINATED  BY  THE  OPERATOR.  IT 
REDUCES  THE  DATA  TO  CALCUALTIONS  OF  NET 
POWER,  RAYLEIGH  NUMBER  AND  NUSELT  NUMBER. 


Emf (70).Pouer(9) . T(60 ) . Tavq( 9) . Ts  (  9 ) 
Tfilm(9).Qnet(9),H(9).K(9).Rho(9).Cp(9) 
N(9).Nu(9).Ra(9).Delt(9).Alfa(9),Pr(9) 
Gr(9) .Beta(9) ,Dpow( 9) , Dts(9) 


COM  /Co/  D(7) 
! 

DIM 
DIM 
DIM 
DIM 
! 

!  CORRLEATION  FACTORS  TO  CONVERT  Emf  TO  REGREES  CELCIUS. 
DATA  0.10086091 ,25727.9.-757345.8,78025596, 
DATA  -9247 486589 , 6 . 98E 1 1 , -2 . b5E 1 3 . 3 . 94E 1 4 
! 

READ  D(*) 

!  RESISTANCE  OF  THE  PRECISION  RESISTOR  IS:  1'/. 

Rp-2.00 

I 

PRINTER  IS  701 

BEEP 

BEEP 

I 

INPUT  -ENTER  THE  NAME  OF  THE  FILE  CONTAINING  DATA' ,OIdf ile$ 

PRINT  USING  "15X."'THE  RAW  Emf  DATA  ARE  FROM  THE  FILE:     "" , lOfi" :01df i le$ 

INPUT  "ENTER  THE  POWER  SETTING  ",Power$ 

PRINT  USING  "17X,""  THE  POWER  SETTING  PER  CHIP  HAS:    "" , 1 OA" ;Power$ 

PRINT 

I 

BEEP 

BEEP 

ASSIGN  ^File  TO  OldfileS 

ENTER  ®File;Emf (») 

I 

!  CONVERT  Emf  TO  DEGREES  CELCIUS 

FOR  1=0  TO  50 

Sum"0 

FOR  J-0  TO  7 

Sun-Sum+D( J)*Emt (I) *J 

NEXT  J 

T(I)»Suni 

NEXT  I 

I 

!  CONVERT  Emf  TO  POWER 

J»1 

Volt»Emf(61) 

FOR  1-62  TO  70 

Pouer(J)-Emf (I)»(Volt-Emf(I))/Rp 

Dpotj(J)-.01*Power(  J) 

J»J+1 

NEXT  I 
I 

!  AREA  OF  THE  BLOCK  FACES 
Acen-1.92E-4 
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530  Alef=1.44E-a 

640  Aria'l .4aE-A 

650  rttoD-^.8E-5 

660  Pibot'M.SE-S 

670  Atot=5.76E-4 

680  ! 

690  ! CALCULATE  THE  AVERAGE  TEMPERATURES  OF  THE  BLOCK  FACES 

700  ! 

710  Tava( I ) • ( T ( 0 ) ♦Acen+T ( 1 ) ♦Atop+TCZ ) ♦Ar la+T ( 3 ) ♦Alef +T ( 4 ) *Abot ) /Atot  - 

720  Tavq(2)-(T(5)*Acen  +  T(7)»Atop-t-T(3)*Arig  +  T(9)*Aief +  T(10)»Abot)/Atot 

730  Tavg(3)-(T( I2)*rtcen  +  T( 1 3 ) ♦Atop  +  T( 1 4 ) ♦Ar la  +  T M 5 ) *A lef *T ( 1 6 ) ♦Abot ) /A  tot 

740  Tavg(4)>(T(  18)*Acen+r(  19)»Atop  +  T(20)*Arig  +  T(21  )  ♦■Alef +  T  (  22  )  ♦Abo  t  )/Atot 

750  Tavql5)=(T(24)»Acen+T(25)»Atop+T(26)*Arig+T(27)*Alef+Tf28)*Abot)/Atot 

760  Tavg(6)  =  (T(30)*Acep*T(31  )  ♦Atop*-T  (  32  )  *-Ar  ig+T  (  33  )  ♦Aief  >T  (  34  )  ^Abot ) /Acq  t 

770  Tava(7)"(T(36)»Acen+T(37)»Atop+T(38>»Arig+T(39)»Alef+T(40)*Abot)/Atot 

780  Tavg(8)  =  (T(42l»Acen  +  T(43)*Atop  +  T(44)*Arig+T(45)*Alef-^T(45)«-Abot  )/Atot 

790  Tavg(9)  =  (T(48)*Acen  +  T(49)*Atop  +  T(50)*Arig  +  T(51  )  ♦Alef  ■••T(52  )  ♦Abot ) /Atot 

800  ! 

810  !  UNCERTAINTY  OF  THE  THERMOCOUPLE  READINGS 

820  Dt=.l 

830  Dtavg>.1 

840  ! 

8S0  !  RESISTANCE  OF  PLEXIGLASS  AND  UNCERTAINTY 

860  Rc=764.S5 

880  ! 

890  !  CHIP  BACK  SURFACE  TEMPERATURES 

900  Ts(1)=T(5) 

910  Ts(2)=T(11) 

920  Ts(3)-T(17) 

930  Ts(4)»T(23) 

940  Ts(5)=T(29) 

950  Ts(5)«T(35) 

960  Ts(7)-T(41) 

970  Ts(8)=T(47) 

980  Ts(9)-T(53) 

990  Tssum"0 

1000  FOR  J=1  TO  9 

1010  Tssum«Tssum+Ts(J) 

1020  NEXT  J 

1030  ! 

1040  Tsava-Tssutn/9 

1050  ! 

1050  !  CONDUCTION  LOSS  CALCULATION  AND  UNCERTAINTY 

1070  Qioss3=(T(17)-Tf55) )/Rc 

1080  Qioss5-(T{29)-T(55) )/Rc 

1090  01oss7-(T(41)-T(54) )/Rc 

1  too  Qloss=(01oss3+aioss5+01oss7)/3 

1110  Dqioss=(Qloss)^.022 

1120  ! 

1130  !  AVERAGE  SINK  TEMPERATURE  CALCULATION 

1140  Tsinlc>(T(57)+T(58)+T(59)*T(60))/4 

1150  Dtsinif-(Dt/T(58)) 

1160  L-8.0E-3 

1170  ! 

1180  PRINT  USING  "IX, '"CHIP      ONET(W)         Tavg-Ts  Nu 

•/.UNC  IN  Nu  '".lOA" 

1130  PRINT 

1200  ! 

1210  ! 

1220  !  CALCULATION  OF  NET  POWER,  Nu,  Ra .  AND  UNCERTAINTIES 
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1230  ! 

1240  FOR  J«1  TO  9 

1241  PRINTER  IS  1 

1242  PRINT  "J" 

1244  PRINT  J 

1245  PRINTER  15  701 
1250  ! 

1260  !  CALULATIQN  OF  Onet 

1270  Onet( J)=Power( J)-Qloss 

1280  Dqnet=(Dpou,lJ)  *2+Dqloss"2)  *  .5 

1290  ! 

1300  !  CALCULATION  OF  Tfilm 

1310  Tf  ilm(J)  =  (Tav/a(J)+Tsinl()/2 

1320  Dtfilm-(Dtsink'2+Dtavq-2)".5 

1330  ! 

1340  !  CALCULATION  OF  A  DELTA  TEMPERATURE 

1350  Delt(  J)  =  Tav/g(J)-Tsinlc 

1360  Ddelt=(Dtavg'2+Dtsink '2) ".5 

1370  ! 

1380  !  CALCULATION  OF  CONVECTION  COEFFICIENT 

1390  H( J)=Onet( J)/(Atot*Delt( J) ) 

1400  Dh=H(J)*( (Dqnet/Onet(J)) *2+(Ddel t/Dei t ( J) )"2)'.5 

1410  ! 

1420  !  CALCULATION  OF  CONDUCTION  COEFFICIENT 

1430  K(J)-(.65-7.89474E-4*Tf  iltn(J))/10 

1440  ! 

1450  !  CALCULATION  OF  DENSITY 

1460  Rho(J)"(1 .825-.00246*Tfilm(J))#1000 

1470  ! 

1480  !  CALCULATION  OF  SPECIFIC  HEAT 

1490  Cp(J)-(  .241in+3.7037E-4»Tf  ilm(J))»4180 

1500  ! 

1510  !  CALCULATION  OF  VISCOSITY 

1520  N(J)-1.4074-2.964E-2»Tfilm(J)+3.8018E-4»Tfilm(J)*2-2.7308E-6*Tfilfn(J)-3+8. 

1679E-9*Tf  ilm(J)-4 

1530  N(J)-N(J)*I .E-S 

1540  ! 

1550  !  CALCULATION  OF  BETA 

1550  Beta(J)-. 00246/(1 .825-.0G246*Tf i lm( J) ) 

1570  ! 

1580  !  CALCULATION  OF  ALPHA 

1590  Alfa(J)-K(J)/(Rho(J)*Cp(J)) 

1600  ! 

1610  !  CALCULATION  OF  PRANDTL  NUMBER  AND  UNCERTAINTY 

1520  Pr(J)=N(J)/Alfa(J) 

1630  ! 

1640  !  CALCULATION  OF  NUSSELT  NUMBER 

1550  Nu(J)-H(J)*L/KCJ) 

1660  Dnu-Nu(J)#(Dh/H(J)) 

1670  ! 

1680  !  CALCULATION  OF  PERCENT  UNCERTAINTY  IN  NUSSELT  NUMBER 

1690  Pernu-(Dnu/Nu(J))»100 

1700  ! 

1710  !  CALCULATION  OF  GRASHGF  NUMBER  AND  UNCERTAINTY 

1720  Gr(J)-9.81*Beta(J)*(L*3)*Deit(J)/N(J)-2 

1730  Dgr-Gr(J)«(Ddelt/Delt(J)) 

1740  ! 

1750  !  CALCULATION  OF  RAYLEIGH  NUMBER  AND  UNCERTAINTY 

1760  Ra(J)-Gr(J)*Pr(J)»1 .E-6 

1770  Dra-Ra(J)*(Dgr/Gr(J))  •■ 

1780  ! 
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179^  i  CALCULATION  OF  FLUX  BASED  RAYLEIGH  NUMBER 

1800  ! 

1310  Raf  (J)  =  '  (9.8I*Beta(  JJ+L'-^+OnetlJ)  )/(  K  (  J  )-»N  (  J  )  *A1  fa  (  J  )  *Atot ))«  1  .  E- 7 

1820  Draf=Raf ( J)»(Dqnet/Qnet(J) ) 

1830  ! 

1840  !  PERCENT  UNCERTAINTY  IN  RAYLEIGH  NUMBERS 

1850  Perra=(Dra/Ral J) J*100 

1860  Perraf=(Drat/Raf( J) )*100 

1870  ! 

1880  PRINT  USING  " 1 X . D . 1 X . 4 ( 7X . ZZ.DDE , ) " ; J . One t ( J ) . Delt ( J) , Nu ( J ) .Pernu 

1890  ! 

1910  PRINT  USING  "6X ," "TEMPERATURE  BASED  RAYLEIGH  NUMBER  ♦  E-G  IS:   '".3D. DDE": 

Ra(  J) 

1920  PRINT  USING  "6X.""7.  UNCERTAINTY  IN  THE  TEMPERATURE  BASED  RAYLEIGH  NUMBER  I 

S  :-".DDD.DDE":Perra 

1330  PRINT  USING  "GX.'-FLUX  BASED  RAYLEIGH  NUMBER  ♦  E-7  IS:   "  "  .  3D . DDE" ;Raf ( J) 

1940  PRINT  USING  "GX.""'/:  UNCERTAINTY  IN  FLUX  BASED  RAYLEIGH  NUMBER  IS:   "",DDD. 

DDE":Perraf 

1950  PRINT 

1960  NEXT  J 

1970  ASSIGN  *File  TO  ♦ 

1980  END 
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3.      FastScan 


\1 

•  PROGRAM  FASTSCAN 

1  1 

20 

!  EDITED  BY  LT.  T  J  BENEDICT 

jn 

! 

40 

!  THIS  PROGRAM  SCANS  3  CHANNELS  FOR  VOLTAGE  MEASUREMENTS 

60 

lDdS9-599 

70 

Pass-0 

71 

N-0 

80 

DIM  T1(S99).V1(2).Y1(S99) 

81 

DIM  T2(S99) .V2(2).Y2(S99) 

82 

DIM  T3(S99).V3(2).Y3(S99) 

90 

CLEAR  709 

100 

CLEAR  122 

101 

1 

102 

BEEP 

103 

PRINTER  IS  701 

lOS 

PEEP 

106 

INPUT  "ENTER  THE  FIRST  FILE  NAME:   ".NewfllelS 

107 

INPUT  'ENTER  THE  SECOND  FILE  NAME:  ".Newfile2$ 

108 

INPUT  -ENTER  THE  THIRD  FILE  NAME:  •.Neufile3$ 

109 

INPUT  "ENTER  THE  VOLTMETER  READING:  '.VS 

113 

PRINT  USING  "ISX, ""THESE  RESULTS  ARE  NOW  STORED  ON  DISK 

11^ 

PRINT 

116 

PRINT  USING  *25X.""FILE:  •" . 1 OA" :Newf i le1$ 

117 

PRINT 

118 

PRINT  USING  ■2SX.""riLE:  ' " . lOA" :Newf i le2$ 

119 

PRINT 

120 

PRINT  USING  •25X.-TILE:  "M  OA"  :Newf  i  le3$ 

121 

PRINT 

123 

WAIT  1 

12^ 

BEEP 

12S 

OUTPUT  709:"AEr 

130 

WAIT  2 

131 

BEEP 

140 

OUTPUT  722:"T4  Fl  R1  PO  ZO  ISTI  SOI  ISTN* 

1A1 

! 

M2 

!  LOOP  NUMBER  ONE 

143 

1 

ISO 

OUTPUT  709: "AFO  ALO" 

160 

OUTPUT  709:-AS" 

170 

BEEP 

17S 

Timedatel-TIMEDATE 

180 

FOR  Jj-0  TO  Ipass 

200 

OUTPUT  722:"T3' 

210 

ENTER  722:V1(») 

230 

TWPass>-Vl  (1) 

250 

Pass-Pass*! 

261 

N-N*1 

260 

NEXT  Jj 

261 

Timedate2-TIMEDATE 

262 

Tota 1 t imel " T imedate2-T imedatel 

263 

OUTPUT  722:'AC0- 

264 

Pass-0 

256 

! 

267 

!  LOOP  NUMBER  TWO 

268 

! 

269 

OUTPUT  709:"AF12  AL12" 

270 

OUTPUT  709:"AS' 

271 

BEEP 

FASTSCAN'    "MOA* 
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272 

274 

27S 

276 

277 

278 

279 

280 

281 

282 

283 

284 

28S. 

286 

287 

288 

289 

291 

292 

293 

294 

296 

296 

297 

298 

299 

300 

otal 

301 

302 

303 

304 

306 

307 

308 

309 

310 

311 

312 

313 

314 

320 

340 

341 

342 
343 
^44 
346 
346 
347 
348 
349 
350 
3S1 
3S2 
363 
364 
366 
366 
357 
358 
390 
400 


PEEP 

rOR  Ji-0  TO  Ipass 

OUTPUT  722;*T3' 

ENTER  722:V2(») 

T2(Pa55) -V2( 1 ) 

Pd55-Pass*l 

NEXT  Jj 

OUTPUT  722:"flC12" 

Pass-0 

LOOP  NUMBER  THREE 

OUTPUT  709:*flF31  AL31" 

OUTPUT  709:'A3* 

BEEP 

BEEP 

BEEP 

FOR  Jt-n  TO  Ipass 

OUTPUT  722:"T3" 

ENTER  722:V3(») 

T3{Pass»-V3(1 ) 

Pass-Pass*1 

NEXT  Jj 


END  LOOPS 


PRINT  USING  ' 
t  imel 
PRINT 

PRINT  USING  ' 
PRINT 

PRINT  USING  ' 
PRINTER  IS  1 
BEEP 


15X,'"THE  TOTAL  TIME  ELAPSED  WAS  (SECONDS) :"  .2X. (DDD. DDE) ' ; 

16X."THE  TOTAL  NUMBER  OF  SCANS  UAS  :   "  "  .2X,  (DDDD.D.2X)  *  ;N 
15X.""THE  VOLTMETER  READING  WAS:   -MOA.'-.VS 


TRANSFER  FIRST  SCAN  DATA 

CREATE  BDAT  Neufile1$.20 

ASSIGN  ftFile  TO  NewfilelS 

OUTPUT  aFile:T1(») 

FOR  Ii-O  TO  Ipass 

T1  (Ii)-. 1008609 1+25727.9*T1(Ii) 

NEXT  Ii 

TRANSFER  SECOND  SCAN  DATA 


CREATE  BDAT  Newfile2$.20 

ASSIGN  «Fi le  TO  Newf ile2$ 

OUTPUT  *Fi le;T2(«) 

FOR  Ii-O  TO  Ipass 

T2(ll)-.1 0086091 ♦25727.9«T2(Il) 

NEXT  1 1 

I 

!  TRANSFER  THIRD  SCAN  DATA 
I 

BDAT  Newfile3$.20 

^File  TO  Newfile3$ 

«File:T3(«) 
■0  TO  Ipass 
■•  10086091 ♦25727. 9*T3(Il)- 767346 


767346.8M1(Ii)-2*78002556«T1(Ii)-3 


767345.8«T2(Il)*2*78002556»T2(Il)'3 


CREATE 
ASSIGN 
OUTPUT 
FOR  Ij- 
T3(ll)' 
NEXT  I 
STOP 
END 


8»T3(Il)*2*78002556«T3(Ii)'3 
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4.   Plot 


i8 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

MO 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

330 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 


FILE  NAME:  PLOT 
EDITED  BY  Lt  T  J  BENEDICT 

THIS  PROGRAM  PLOTS  THE  DATA  ACQUIRED  BY 
PROGRAM  "FASTSCAN". 

PRINTER  IS  705 

BEEP 

Xmin-0 

Xmax-200 

BEEP 

INPUT  "ENTER  MINIMUM  AND  MAXIMUM  Y-VALUES"  ,  Ymin. Ymax 

BEEP 

Xstep-20 

BEEP 

Ystep-.2 

BEEP 

PRINT  '•IN:SP1:IP  2000.2000.8000,7000:" 

PRINT  'SC  0. 100,0. 100:TL  2,0;' 

St  x-100/(Xmax-Xmin) 

SfylOO/(Ymax-Ymin) 

PRINT  "PU  0.0  PD" 

FOR  Xa-Xmin  TO  Xmax  STEP  Xstep 

X-(Xa-Xmin)*Sfx 

PRINT  "PA";X,',0;  XT;" 

NEXT  Xa 

PRINT  "PA  100.0:PU:" 

PRINT  "PU  PA  0.0  PD' 

FOR  Ya-Ynin  TO  Ymax  STEP  Ystep 

Y-(Ya-Ymin)*Sfy 

PRINT  "PA  0,";Y,"YT" 

NEXT  Ya 

PRINT  "PA  0,100  TL  0  2" 

FOR  Xa-Xnin  TO  Xmax  STEP  Xstep 

X'(Xa-Xmin)*Sfx 

PRINT  •PA';X,",100;  XT" 

NEXT  Xa 

PRINT  "PA  100,100  PU  PA  100,0  PD" 

FOR  Ya-Ymin  TO  Ymax  STEP  Ystep 

Y»(Ya-Ymin)»Sty 

PRINT  "PD  PA  100,",Y,"YT" 

NEXT  Ya 

PRINT  "PA  100.100  PU" 

PRINT  "PA  0,-2  SR  1 .5,2" 

FOR  Xa=Xmin  TO  Xmax  STEP  Xstep 

X" ( Xa-Xmin) ♦Sf x 

PRINT  "PA";X.,",0;" 

PRINT  "CP  -2.-i;LB":Xa;"" 

NEXT  Xa 

PRINT  "PU  PA  0.0" 

FOR  Ya-Ymin  TO  Ymax  STEP  Ystep 

IF  ABS(Ya)<1 .E-5  THEN  Ya-0 

Y-(Ya-Ymin)*Sfy 

PRINT  "PA  0.";Y."" 

PRINT  "CP  -5,-.25;LB";Ya;"" 

NEXT  Ya 

BEEP 

IdI-0 

IF  Idi-0  THEN 
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610  BEED 

520  KiaDel5="Time  (sec)" 

630  BEEP 

640  Ylabel$="TempGrature  (C)' 

650  PRINT  "SR  1.5.2:PU  PA  50.-10  CP" ; -LENf XlabelS  ) /2; "0 :LB- : XlabelS : 

650  PRINT  "PA  -11.50  CP  0 . " ; -LEN( Y labe  IS ) /2»5/5  ;  "DI  0 . 1 :LB" : Y labels; 

670  END  IF 

680  PRINT  "CP  0,0' 

690  BEEP 

700  INPUT  "ENTER  THE  NAME  OF  THE  DATA  FILE" .D_f i leS 

710  ASSIGN  *File  TO  D_fiie$ 

720  BEEP 

730  Md=0 

740  BEEP 

750  Npairs=500 

760  BEEP 

770  PRINTER  IS  1 

780  Sym'1 

790  PRINTER  IS  705 

800  PRINT  "PU  DI" 

810  IF  Sym=1  THEN  PRINT  "SM." 

820  IF  Sym=2  THEN  PRINT  •SM+" 

830  IF  Sym=3  THEN  PRINT  "SMo" 

840  IF  Md>1  THEN 

850  FOR  1-1  TO  (Md-1) 

860  ENTER  *File;Xa.Ya 

870  NEXT  I 

880  END  IF 

890  FOR  Xa-0  TO  199  STEP  .3333333 

900  ENTER  ^File;Ya 

910  Ya-.10085091>25727.9»Ya-7S73^5.8»Ya'2+78002556*Ya'3 

920  X»(Xa-XmiTi)*Sfx 

930  Y-(Ya-Ymin)*Sfy 

940  IF  Sytn>3  THEN  PRINT  "SM" 

950  IF  Sym<4  THEN  PRINT  "SR  1.4,2.4" 

960  PRINT  'PA'.X.Y.'PD* 

970  IF  Sym>3  THEN  PRINT  "SR  1.2,1.6" 

980  IF  Sym=4  THEN  PRINT  •UC2, 4 ,99, 0 , -8. -4 , 0 , 0 .8, 4 . 0 . ; " 

990  IF  Sym=5  THEN  PRINT  "003.0 ,99, -3, -6, -3.6.3.6.3 , -6; " 

1000  IF  Sym-5  THEN  PRINT  "UCO .5. 3.99. 3. -8 . -6 . 0 .3.8: " 

1010  IF  Sym=7  THEN  PRINT  "UCO ,-5 .3 ,99, -3 ,8 .6, 0 , -3, -8: " 

1020  NEXT  Xa 

T021  PRINT  "PU" 

1030  BEEP 

1040  ASSIGN  ?>FiIe  TO  ♦ 

1050  END 
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